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ABSTRACT 


The  temperature  of  formation  of  the  gold-quartz  deposits,  Yellow¬ 
knife,  N.W.T.,  is  obtained:  1)  by  studying  isotope  fractionation  among 
the  coexisting  sulfide  minerals  associated  with  gold-quartz  veins,  and 
2)  by  fluid  inclusion  homogenization  temperatures  on  quartz,  to  be  in 
the  range  of  350°-300°C. 

From  the  thermochemical  data  and  the  knowledge  of  the  mineral 

assemblages,  the  fugacities  of  oxygen  and  sulfur  are  calculated  to  be 

in  the  range  of  10  exp.  -31  to  -29  and  -10.5  to  -8.5  respectively. 

18 

The  0  of  the  hydrothermal  fluids  is  calculated  assuming  equi- 

18 

librium  between  quartz-calcite-water  (0  =  5.0-6.8%o)  and  this 

water 

value  is  very  close  to  the  established  isotopic  compositions  of  magmatic 
waters . 

From  the  derived  fugacities  of  CO^  and  CH^ ,  the  total  carbon 
isotope  values  are  calculated  to  be  in  the  range  of  -4.2  to  -0.5%o, 
based  on  the  mole  fractions  of  CO^  and  CH^ ,  assuming  CaC) ^-CO^-CH^ 
equilibria. 

By  comparing  the  frequency  plot  of  Na/K  ratios  obtained  on  the 
fluid  inclusions  from  quartz  with  the  established  plots  of  different 
aqueous  solutions  the  ore-forming  fluids  are  concluded  to  have  been 
derived,  most  likely,  from  a  magmatic  source. 
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INTRODUCTION 

Yellowknife  is  situated  on  the  northern  shore  of  Great  Slave 
Lake  (lat  62°30’N,  long  114°30TW)  in  the  North  West  Territories.  The 
Con-Rycon  Mine  is  located  nearly  2  miles  South  of  the  Yellowknife 
town  site  (Fig.  1).  The  Con-Rycon  Mine,  which  began  production  in 
1938,  is  the  major  gold  producer  of  the  district  and  is  the  property 
of  Consolidated  Mining  and  Smelting  Company  of  Canada.  The  other 
producing  mine  operating  in  the  area  is  the  Giant  Mine,  situated 
about  3  miles  North-east  of  the  Con-Rycon  Mine. 

The  topography  of  the  Yellowknife  area  is  typical  of  the  Pre- 
Cambrian  Shield.  The  country  is  almost  flat.  Numerous  lakes  and 
muskegs  are  quite  common  in  the  area.  Permafrost  conditions  exist 
in  most  parts  of  the  area. 

Objectives  of  the  Thesis 

Laboratory  investigations,  including  optical,  isotopic  and  fluid 
inclusion  studies  have  been  conducted  on  the  samples  collected  from  the 
Campbell  Shear  System,  Yellowknife.  It  is  the  purpose  of  the  thesis  to 
reconstruct  a  possible  geological  environment  in  the  light  of  the 
physico-chemical  factors  that  would  have  influenced  ore  deposition. 

The  techniques  used  in  this  study  were:  (a)  sulfur  isotope  studies, 
(b)  fluid  inclusion  work,  (c)  carbon  and  oxygen  isotope  work,  (d)  K/Na 
ratio  distributions  in  fluid  inclusions,  (e)  mineral  associations,  and 
(f)  application  of  the  principles  of  thermodynamics  in  natural  systems 


and  isotope  geochemistry. 
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CHAPTER  I 


REGIONAL  GEOLOGY 


The  mineral  potential  of  the  Yellowknife  area  has  resulted  in 
extensive  detailed  mapping  of  this  district  of  the  Slave  Province  since 
1938  (Henderson,  1939;  Henderson  and  Brown,  1949,  1952,  1966;  Lord, 

1942;  Fortier,  1947;  Folinsbee,  1947;  Moore,  1956;  Boyle,  1961;  Baragar, 
1966;  and  Kretz,  1967). 

The  Precambrian  rocks  of  the  Yellowknife  group  have  been  broadly 
divided  into  two  major  rock  units:  Division  A  and  B  (Joliffe  1938). 
Division  A  consists  of  a  series  of  dark  green,  metamorphosed  Volcanic 
rocks  which  includes  meta-andesites,  meta-basalts  and  intercalated  meta- 
dacites  with  associated  pyroclastics .  Division  B  consists  predominantly 
of  conglomerates,  arkosic  quartzites  (sub-greywackes)  and  acid  volcanic 
rocks  resting  more  or  less  unconformably  on  the  rocks  of  Division  A 
(Table  1) .  Palimpsest  volcanic  structures  are  often  preserved  in  the 
rocks  of  Division  A,  especially  pillow  and  flow-structures,  which  are 
quite  noticeable  on  the  surface  as  well  as  at  the  underground  levels  of 
the  mines.  Slight  modifications  in  categorizing  the  various  rocks  have 
been  introduced  by  later  authors  (Henderson  and  Brown,  1966;  Green,  1968), 
however,  the  major  distinction  between  Divisions  A  and  B  has  been 
retained  in  all  cases.  The  volcanic  and  sedimentary  rocks  (Division  A 
and  Division  B)  are  intruded  by  many  sills  and  dykes  of  meta-gabbroic 


and  meta-dioritic  character. 


FIGURE  1 


Regional  Geology,  Yellowknife  Area, 
Northwest  Territories 
(After  Boyle  1961) 


MAP  ADAPTED  FROM  BOYLE,  1961,  FlG.I. 
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TABLE  1 


Table  of  Formations  (Pre-Cambrian) 
Yellowknife  Area,  Northwest  Territories 
(adopted  from  Green  1968) 


FORMATION 


DESCRIPTION 


REMARKS 


AGE 


Diabase  dykes 

Granite,  granodiorite 
Nepheline-sodalite  syenite 

(Hudsonian) 

PROTEROZOIC 

Diabase,  picrite-gabbro 
dykes 

(Leech  1966) 

Snare  Group 

Finely  bedded,  alternating 
siltstone  and  shale  and 
metamorphic  equivalents. 
Dolomite,  arkose 

(Lord  1942) 
(Ross  and 
McGlynn  1963) 

UNCONFIRMITY 

Redout  Lake,  Prosperous 

Lake  granites  and  associ¬ 
ated  pegmatites 

Western  granodiorite 
(?  Ross  Lake  granodiorite) 
Stock  Lake  stock,  aplites 
(?  South-east  granodiorite) 

ARCHAEAN 

Quartz  feldspar  porphyries 

Second  phase  of  metadiorite 
and  meta  gabbro  intrusion 

Yellowknife  Group 
Division  B 

Graywacke,  argillite, 
phyllite,  quartzite  etc. 
and  their  metamorphic  equi¬ 
valents.  Dacite,  trachyte 
rhyolite,  quartz  porphyry, 
pyroclastics  and  their 
metamorphic  equivalents. 

(Overlie 

Div.  A  with 
apparent 
regional 
conf irmity) 

Unclassified 
(stratigraphic  position 
uncertain) 

Conglomerate  and  quartzite 

Local  Unconformity 

?  (South-east  granodiorite)  (Folinsbee  et  al.  1968) 

First  phase  of 
and  metagabbro 

metadiorite 

intrusion 

Yellowknife  Group 
Division  A 

Basalt,  andesite,  dacite 
rhyolite,  cherty  tuff 
and  their  metamorphic 
equivalents 

(Ross  Lake  Granodiorite)  (Baragar  1966) 

1.1  METAMORPHIC  FACIES  -  YELLOWKNIFE  GROUP 


The  volcanic  rocks  of  the  Yellowknife  group  have  been  considerably 
altered  during  regional  metamorphism.  These  basic  volcanics  exhibit 
three  metamorphic  grades,  the  isograds,  which  are  almost  parallel  to 
the  Western  granodiorite  greenstone  contact  (see  Fig.  1) .  The  amphi¬ 
bolite  facies  rocks  adjoining  the  Western  granodiorite  grade  into 
epidote-amphibolite  lithologies  followed  by  rocks  of  green  schist  facies, 
when  one  traverses  from  West  to  East. 

The  rocks  of  amphibolite  facies  grade  are  mainly  fine-grained, 
meta-basalts  and  meta-andesites.  Horenblende,  plagioclase  (andesine) 
and  epidote  are  the  major  constituent  minerals.  Quartz,  ilmenite  and 
sphene  occur  in  small  quantities.  In  contrast  the  rocks  of  the  epidote 
amphibolite  facies  are  relatively  massive  and  contain  light-green, 
fibrous  amphibole,  plagioclase  (albite  to  oligoclase)  and  epidote. 

Quartz  and  small  amounts  of  carbonate  minerals  are  also  almost  ubiquitous. 

The  rocks  of  the  green  schist  facies  are  of  relatively  limited 
occurrence  in  the  district.  One  of  the  areas  where  these  rocks  are 
exposed  is  the  so-called  Giant-fault  block  East  of  Bow  Lake.  The 
constituents  of  the  green  schist  lithologies  are  green  chlorite,  epidote 
plagioclase  (albite)  and  carbonate  minerals.  The  most  characteristic 
feature  of  the  meta-volcanics  is  that  they  have  often  retained  the 
original  pillow  structures.  According  to  Henderson  and  Brown  (1966) 
there  are  two  varieties  of  volcanics:  massive  flows  and  pillowed  flows. 
Gradations  between  these  two  varieties  are  often  seen.  The  massive  flows 
are  relatively  coarser  grained.  The  dimensions  of  the  pillows  range  in 
length  from  a  few  inches  to  more  than  1-2  feet.  Each  pillow  is  separated 
from  its  neighbouring  pillow  by  a  strong  alteration  rim  surrounding  each 
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pillow. 

Henderson  and  Brown's  (1966)  petrographic  study  and  Boyle's 
(1961)  chemical  studies  on  these  meta-volcanics  may  be  summarized  as 
follows:  SiO^  and  CaO  show  little  difference  from  a  typical  basalt, 
whereas  and  show  a  slight  depletion  and  H^O,  CO^,  K^O,  and 

Na^O  exhibit  a  minor  enrichment  above  a  normal  basaltic  composition. 
There  are  different  opinions  on  the  origin  of  the  pillowed  structure  in 
the  pre-Cambrian  volcanics.  None  of  the  theories  would  provide  a  clear 
explanation  of  the  origin  of  the  pillow  structures.  Lewis  (1941) 
suggested  that  these  pillow-lavas  have  originated  due  to  sub-aqueous 
conditions.  Recent  views  of  Henderson  and  Brown  (1966)  on  the  mode  of 
formation  of  the  pillows  are  that  each  pillow  has  been  emplaced  as  an 
individual  unit,  and  in  most  cases,  the  upper  pillow  is  deposited  when 
the  lower  pillow  is  still  in  semi-solid,  flexible  plastic  state.  It 
is  believed  by  most  of  the  authors  that  pillow  structure  would  probably 
have  resulted  from  sudden  chilling  action  at  the  contact. 

Meta-dacites  and  associated  meta-agglomerates,  meta-tuffs  and  meta¬ 
breccia  are  exposed  West  of  the  old  town-site,  near  the  West-Bay  fault. 
Apparently  these  flows  are  displaced  North  to  South  by  more  than  three 
miles  by  the  West-Bay  fault  and  their  continuation  could  be  seen  at  the 
Giant  C  shaft.  These  dacites  are  generally  fine-grained,  white  to  buff- 
colored  when  weathered  and  often  contain  phenocrysts  of  feldspar  and 
quartz.  Henderson  and  Brown  (1966)  have  pointed  out  that  it  is  often 
difficult  to  assess  the  exact,  original  mineralogical  composition  of 
these  rocks.  In  Henderson's  opinion  (1966),  the  more  siliceous  and 
porphyritic  types  might  have  been  derived  from  original  normal  basalts 
and  the  meta-andesites  could  be  the  result  of  later  processes  of  alter- 
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ation  and  metasomatism  involving  albitisation ,  silicif ication  and 
chloritization. 

1.2  PLUTONIC  INTRUSIVES 

Certain  plutonic  intrusive  masses  in  the  district  have  been  studied 
in  detail  by  Green  (1968) .  The  Western  granodiorite  is  a  relatively 
large  intrusive  body  which  abuts  against  the  greenstone  belt  and  extends 
more  than  50  miles  westwards.  The  composition  of  this  intrusive  varies 
from  granite,  adamellite,  granodiorite  to  quartz  diorite.  Compositional 
variations  in  Na/K  ratios  have  been  recorded  by  Boyle  (1961)  who  inter¬ 
preted  the  higher  Na/K  ratios  along  the  eastern  margin  of  the  intrusive 
body  (compared  with  the  ratios  in  the  interior  portions)  as  being  due  to 
K+  mobilisation  into  the  metavolcanics .  In  the  northern  part  of  the 
area,  the  contact  between  the  granodiorite  and  metavolcanics  is  grad¬ 
ational  whereas  in  the  southern  parts  it  is  quite  sharp  and  irregular. 
Green’s  modal  analyses  (1968)  shows  a  gradational  variation  from  granite 
to  granodiorite  and  adamellite.  The  Stock  Lake  granodiorite  is  similar 
to  the  Western  granodiorite  in  texture  and  composition  (Henderson  and 
Brown,  1966).  There  is  clear  evidence  of  assimilation  of  the  bordering 
greenstones  by  the  intrusives  at  their  contacts.  The  composition  of 
the  Stock  Lake  granodiorite  varies  from  quartz  diorite  to  quartz  monzonite 
(Green  1968)  .  The  southeast  granodiorite  ranges  in  composition  from 
normal  granodiorite  to  trondhjemite  (Green  1968)  but  it  is  predominantly 
biotite-bearing  granodiorite.  The  Ross  Lake  granodiorite,  the  Redout 
Lake  granite  and  the  Prosperous  Lake  granite  are  the  other  intrusive 
bodies  identified  in  the  area,  and  their  relationships  with  the  metavolcanics 


remain  to  be  worked  out  in  detail. 
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CHAPTER  II 

GEOLOGY  OF  THE  SHEAR  ZONES  AND  ORE  DEPOSITS 

Detailed  descriptions  of  the  shear  zones  and  the  nature  of  the  ore 
deposits  have  been  provided  by  Boyle  (1961)  and  this  work  serves  as  the 
basis  for  the  descriptions  and  discussions  embodied  in  this  thesis. 

The  shear  zones  in  the  greenstone  belt  apparently  fall  into  two  structural 
groups,  namely  (a)  those  which  parallel  the  strike  of  the  volcanic  flows 
and  tuff  beds,  and  (b)  those  which  transect  these  horizons.  The  shear 
zones  paralleling  the  flows  and  tuffs  are  narrow,  discontinuous  entities 
which  do  not  carry  any  economic  deposits.  In  contrast,  those  shear 
zones  which  transect  the  greenstone  belt  often  contain  potential  economic 
ore  deposits. 

Boyle's  detailed  descriptions  (1961)  of  the  shear  zones  and  their 
internal  nature  may  be  summarised  as  follows:  two  types  of  internal 
features  may  be  distinguished:  (a)  Breccia  zones,  and  (b)  Schist  zones. 
Breccia  zones  which  are  well  developed  in  the  Negus-Rycon  and  con 
systems  are  characterised  by  lithologies  broken  into  spindle-shaped 
fragments  locally  cemented  by  carbonates.  The  fragments  are  rarely 
greater  than  6  inches  in  any  dimension  and  the  majority  of  them  are  about 
the  size  of  brazil  nuts,  which  they  simulate  in  shape.  The  fragments 
in  these  breccias  are  massive  to  slightly  schistose  in  character  and 
seldom  exhibit  the  marked  schistosity  or  extreme  alteration  of  the  rocks 
which  comprise  the  schist  zones.  Where  quartz  lenses  occur  in  this  type 
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of  shear  zone,  a  buff  colored,  dense  and  dyke-like  alteration  zone  is 
present  in  which  the  outlines  of  the  original  breccia  fragments  can  be 
discerned  upon  close  examination. 

The  schist  zones  are  characterized  by  linear  zones  of  extremely 
schistose  rock  with  very  strong  lineations,  characterized  by  abundant 
chlorite  and/or  sericite.  Schistose  zones  are  especially  common  in  the 
ore  bodies  of  the  Con  and  Negus-Rycon  systems  and  are  the  principal 
features  of  the  Giant-Campbell  shear  zone  system.  In  nearly  all  cases 
the  dips  and  in  some  cases  the  strikes  of  the  schistose  planes  intersect 
the  walls  of  the  shear  zones  at  a  small  angle. 

In  most  of  the  transecting  shear  systems,  gradational  contacts 
between  schist  zones  and  breccia  zones  are  common.  This  feature  is  well 
illustrated  in  the  Negus-Rycon  system  where  the  gradation  can  be  traced 
with  some  degree  of  confidence.  In  some  of  the  shear  zones,  the  ideal 
breccia  shear  zones  grade  imperceptibly  to  schist  zones  along  strike 
and  in  some  cases  down  dip.  In  the  transition  from  the  breccia  type 
zone  to  schist  zones  the  fragments  of  the  breccia  become  elongated  and 
flattened,  alteration  increases  concomitantly  and  the  fractures  separating 
the  fragments  are  partly  or  completely  obliterated.  The  end  result  is  a 
well-developed  schist  in  which  traces  of  fragments  may  appear,  but  in 
general  no  hint  of  the  original  breccia  remains.  The  transition  from 
breccia  to  schist  can  be  seen  in  schist  zones,  but  in  others  no  such 
genetic  origin  for  the  schist  zone  can  be  evoked.  The  internal  features 
of  the  schist  zones  paralleling  the  flows  and  those  of  the  Giant-Campbell 
and  portions  of  the  Con  systems  suggests  that  extensive  brecciation  was 
not  a  major  factor  in  their  development. 
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2.1  THE  CON-RYCON  SYSTEM 

Detailed  descriptions  of  the  main  shear  zones  have  been  supplied 
by  Campbell  (1949) ,  Lord  (1951)  and  Sproule  (1952) .  The  structure  has 
been  traced  southwest  along  the  southeast  end  of  Pud  Lake  and  up  to  the 
Kam  fault  zone  (Fig.  1).  Brown  and  Henderson  stated  that  the  Con 
system  has  branched  out  into  two  or  more  zones  and  is  a  part  of  the 
Campbell  system.  However,  no  positive  evidence  for  such  a  suggestion 
is  apparently  available  at  present  and  any  conclusions  remain  speculative. 
On  the  other  side  of  the  Pud  Lake,  the  Con  shear  system  exhibits  quite 
distinct  characteristics.  On  the  southern  part  of  the  shear  zone,  it  is 
characterized  by  the  coalescence  of  numerous  shear  zones  whereas  in  the 
northern  part  of  the  system,  it  exhibits  a  multiple  branching  effect 
into  series  of  small  shear  zones.  But  on  the  500  ft  level  of  the  Con 
mine  all  these  structures  merge  to  form  the  so-called  C-4  schist  zone 
and  again  at  deeper  levels  the  shear  system  splits  into  minor  shear  zones 
separated  by  several  minor  blocks  of  barren  rneta  andesitic  lithology.  As 
Henderson  and  Brown  (1966)  pointed  out,  the  ore  bodies  are  quartz  lenses, 
pods  and  replacement  bodies  mineralized  with  pyrite,  arsenopyrite ,  chalco- 
pyrite,  sulfo  salts,  galena  and  native  gold. 

The  general  strike  of  the  system  is  N80°E  and  the  dip  53°W.  The 
Negus-Rycon  system  strikes  N10-250  and  dips  45-65°W.  Henderson  and 
Brown  (1966)  suggest  that  the  Negus-Rycon  system  was  probably  formed  by 
thrust  movements  generated  by  the  same  forces  which  produced  the  Giant 
Campbell  and  Con  shear  systems.  The  veins  are  quite  discontinuous,  small 
and  contain  a  wide  variety  of  metallic  minerals:  e.g.,  pyrite,  gold, 
arsenopyrite,  chalcopyrite ,  stibnite,  various  sulfo  salts,  sphalerite  and 


galena. 
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2.2  THE  GIANT -CAMPBELL  SYSTEM 

Extensive  work  has  been  performed  by  many  authors  on  different 
aspects  of  the  Campbell  shear  system  (e.g.  Henderson  and  Brown,  1952b; 
Dadson  and  Bateman,  1948;  Dadson,  1949;  Campbell,  1949;  Lord,  1951; 
Bateman,  1952;  Boyle  et  al.  1961;  and  Henderson  and  Brown,  1966).  The 
Campbell  shear  system  is  one  of  the  district’s  longest  ore-producing 
systems.  The  northern  extensions  of  the  Campbell  shear  system,  on  the 
eastern  side  of  the  West  Bay  fault,  have  posed  a  controversial  problem 
to  many  geologists.  Campbell’s  (1949)  location  of  the  extensions  of 
the  shear  zone  on  the  southern  part  of  the  West  Bay  fault  was  a  note¬ 
worthy  step  to  economic  development  of  the  area. 

The  southern  extension  of  the  Giant-Campbell  system,  south  of 
the  West  Bay  fault  is  known  as  the  Campbell  shear  system  and  is  exposed 
for  8000  ft  in  the  underground  workings  of  the  Negus  and  Con  mines  at 
the  2300  ft  level.  Examination  of  the  mine  plans  and  sections  at 
varous  levels,  indicates  that  the  width  of  the  shear  zone  at  the  2300  ft 
level  is  in  excess  of  200  ft,  whereas  at  the  4500  ft  level  the  width  is 
only  around  95  ft  (Fig.  2).  Moreover,  the  shear  zone  at  the  2900  ft  level 
is  quite  complex  with  intertwined  zones  separated  by  a  barren  zone  of 
massive  andesite.  The  massive  included  block  of  andesite  is  almond 
shaped.  The  shear  zone  maintains  a  strike  direction  of  N25°E  and  dips 
vary  from  hanging  wall  to  foot  wall.  The  foot  wall  dips  are  generally 
steeper  and  the  hanging  wall  dips  vary  from  45°  to  65°W. 

In  the  Campbell  shear  system,  drag-folded  and  contorted  schist 
zones  apparently  acted  as  susceptible  locations  for  accumulating  hydro- 
thermal  solutions  which  induced  silicif ication ,  sericitization  and 


metallic  mineralization. 


FIGURE  2 


Typical  Geological  Section  of  the  Campbell  Shear  Zone 
showing  the  relationship  between  the  Gold-Quartz  Lenses 

and  the  Alteration  Zones 


—  .  . 
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The  mineralogy  of  the  schist  zone  is  comparatively  complex,  being 
composed  of  quartz,  albite,  carbonate  and  sericite. 

2.3  QUARTZ  LENSES  AND  THE  VEIN  PATTERN  IN  THE  SHEAR  ZONE 

Quartz  lenses  enclosed  within  an  envelope  of  sericite-chlorite 
schist  constitute  the  main  ore  bodies  in  the  shear  zone.  Two  generations 
of  quartz  have  been  identified  (Boyle  1961) .  The  first  generation  of 
quartz  occurs  as  complex,  gold-bearing  lenses  and  veins.  The  second 
type  occurs  as  small  stringers  cutting  the  shear  zone.  The  gold-bearing 
lenses  range  from  a  few  inches  in  width  to  3  ft  or  more.  However, 
smaller  en  echelon  lenses  are  more  common  and  these  only  occasionally 
merge  to  form  larger,  richer  lenses.  Small  stringers  of  carbonates 
cutting  across  the  quartz  lenses  as  well  as  the  schistosity  of  the 
chlorite  schist  and  sericite  schist  are  quite  common.  The  nature  and 
mineralogy  of  the  alteration  zone  are  discussed  elsewhere. 

Many  different  opinions  exist  concerning  the  origin  of  the  shear 
zone.  The  internal  structures  of  the  schist  zones,  particularly  the  ore- 
bearing  sericite  and  quartz  zones,  show  clear  indications  of  folding 
(Brown,  Dadson  and  Wriggleworth ,  1959)  which  might  be  the  result  of 
simple  folding  of  the  green  schist  belt.  Secondly  the  schist  zone  might 
have  originated  due  to  shearing  rather  than  folding,  as  the  schist  zone 
transects  the  regional  strike  of  the  flows  at  a  low  angle. 

The  Giant-Campbell  system,  a  complex,  intricate  system  of  shear 
zones  (Henderson  and  Brown,  1966)  is  visualized  as  a  result  of  the 
principal  folding  process.  It  is  believed  that  ore  deposition  accompanied 
the  last  stages  of  development  of  the  shear  zone  system.  No  definite 
metamorphic  or  structural  event  for  the  development  of  dilated  zones  in 
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the  shear  zone  has  been  pointed  out,  however,  it  is  suggested  that  the 
shear  zone  has  originated  as  a  result  of  thrust  faulting  along  which 
the  shear  zone  has  developed  (Henderson  and  Brown,  1966) . 
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CHAPTER  III 

MINERALOGY  AND  PARAGENESIS 

The  principal  objective  of  the  microscopic  examination  of  the 
specimens  discussed  herein  was  the  study  of  textural  relationships 
between  the  co-existing  phases  in  various  samples  collected  at  different 
elevations  in  the  mine,  thereby  interpreting  the  stages  of  deposition, 
the  paragenetic  sequence  and  the  physico-chemical  conditions  prevailing 
during  the  period  of  mineral  deposition.  At  the  same  time  the  micro¬ 
scopic  examination  of  the  ores  served  to  check  the  purity  of  the  mineral 
separates  utilized  in  isotope  analyses  and  to  indicate  the  optimum  size 
for  the  separation  of  pure  fractions.  X-ray  methods  were  concluded  to 
be  inapplicable  to  these  samples  owing  to  the  very  low  percentages  of 
contamination,  likely  to  be  encountered  in  the  mineral  separates  used 
for  analysis. 

The  ore  minerals  identified  were  in  order  of  abundance:  pyrite, 
galena,  sphalerite,  arsenopyrite ,  chalcopyrite ,  native  gold  and  pyrrhotite. 
The  major  gangue  minerals  were  found  to  be  quartz,  calcite,  sericite  and 
chlorite.  Especially  in  the  mineralised  lenses,  quartz  is  the  major 
phase,  sulfides  representing  5  to  8%  by  volume.  Among  the  sulfides, 
pyrite  constitutes  70-80%  of  the  fraction. 
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3.1  DESCRIPTION  OF  THE  MINERALS 
Pyrite: 

Pyrite  occurs  as  large  euhedral  to  subhedral  grains  ranging  from 
0.5  to  2.0  mm  in  diameter.  Most  of  the  specimens  contain  euhedral  pyrite 
crystals  exhibiting  no  replacement  textures.  However,  the  pyrite,  in 
some  instances  does  exhibit  replacement  textures  in  which  pyrite  is 
replaced  by  galena,  sphalerite  and  chalcopyrite .  Arsenopyrite  does  not 
seem  to  have  replaced  pyrite  in  any  specimens. 

Arsenopyrite : 

Arsenopyrite  occurs  mainly  as  euhedral  grains.  It  is  generally 
believed  by  many  of  the  mine  geologists  that  wherever  arsenopyrite  is 
present,  higher  values  of  gold  should  be  expected.  Actually  there  does 
not  seem  to  be  any  genetic  relationship  between  gold  and  arsenopyrite 
and  in  a  mumber  of  instances,  native  gold  was  identified  associated  with 
quartz  and  pyrite  in  sections  devoid  of  arsenopyrite.  The  general  grain 
size  of  arsenopyrite  ranges  from  0.2  mm  to  1  mm  diameter. 

Galena : 

Galena  was  always  associated  with  sphalerite  and  pyrite  and  in 
parts  the  mineral  was  observed  replacing  pyrite.  However,  other  specimens 
showed  alternating  bands  of  galena  and  pyrite  in  apparent  equilibrium. 

The  galena  grains  were  euhedral  to  subhedral,  0.5  mm  to  2  mm  in  diameter. 

Chalcopyrite: 

Chalcopyrite  is  in  most  cases  associated  with  pyrite,  essentially 
replacing  pyrite.  Generally  the  grain  size  of  the  chalcopyrite  ranges 


from  0.2  to  1  mm  and  the  mineral  is  subhedral. 
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Pyrrhotlte: 

Pyrrhotite  was  also  noted  in  a  few  sections  occurring  as  grains 
always  associated  with  pyrite.  It  is  possible  that  pyrrhotite  may  be 
cogenetic  with  the  pyrite. 

Gold: 

Gold  is  apparently  present  in  two  distinct  phases:  (a)  in  the  form 
of  the  native  element,  and  (b)  in  the  form  of  auro  -complexes.  Native 
gold  could  be  identified  very  easily  in  hand  specimens.  Some  of  the 
alteration  zones  as  well  as  the  lenses  show  very  high  assay  values,  even 
though  no  visible  gold  is  present.  In  these  alteration  zones,  most  of 
the  gold  is  probably  present  in  sulf o-complexes . 

Quartz : 

Boyle  (1961)  distinguished  two  types  of  quartz  in  the  ore  bodies 
of  Yellowknife  region;  viz.  earlier  quartz  which  is  massive,  white  in 
colour  and  commonly  occurs  in  small,  narrow,  irregular  lens-like  ore 
bodies  grading  into  the  alteration  zones.  A  later  quartz  occurs  in  small 
and  narrow  veinlets  or  stringers  cutting  earlier  quartz  and  alteration 
zones . 

In  thin  sections  the  quartz  is  medium  to  coarse  grained  and  most  of 
the  sections  show  highly  strained  and  crushed  grains  with  intermittent 
clear  patches  of  unstrained  quartz  between  them.  The  fluid  inclusions 
in  quartz  are  mainly  of  two  sizes,  the  majority  being  0.05  mm  to  0.1  mm 
in  diameter  whereas  others  range  from  0.2  to  0.3  mm.  The  smaller 
inclusions  present  considerable  difficulties  to  determine  the  homogenization 


temperatures . 
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Calcite: 

Calcite  is  the  major  carbonate  mineral  occurring  in  the  quartz- 
carbonate  lenses  as  well  as  the  sericite-carbonate  schist  of  the  alter¬ 
ation  zones.  Many  phases  of  carbonates  have  been  identified  (Boyle  1961). 
Other  carbonates  present  are  ankerite  and  dolomite:  the  fluid  inclusions 
in  these  carbonates  are  very  small  and  it  is  most  difficult  to  obtain  a 
consistant  temperature  determination  with  calcite.  The  calcite  exhibits 
a  well  crystalline  nature  and  is  transparent. 

3.2  PARAGENESIS 

A  description  of  the  ores  from  different  shear  zones  was  given  by 
Ridland  (1941)  and  by  Graham  and  Kaiman  (1951) .  These  authors  proposed 
that  the  mineralization  within  the  shear  zone  might  have  taken  place  in 
many  distinct  stages.  Ridland  identified  33  mineral  species  and  in  his 
opinion  the  mineralogy  of  the  ores  is  quite  complex.  The  mineral  species 
identified  by  Ridland  (1941)  were  gold,  pyrite,  sphalerite,  galena,  chalco- 
pyrite,  arsenopyrite ,  chalcocite,  covellite,  gudmundite,  tenantite, 
pyrrhotite,  leucopyrite,  cobaltite,  stibnite,  nagyagite,  sylvanite,  tan- 
talite,  altaite,  calaverite,  rickardite,  malachite,  siderite,  limonite 
and  hematite.  Sulfosalts  recognised  were  boulangerite ,  meneghinite, 
jamesonite,  berthierite,  bournonite  and  tetrahedrite.  Coleman  (1957) 
postulated  that  gold  mineralization  took  place  during  the  latter  stages 
of  mineralization.  As  can  be  seen  from  Fig.  3,  pyrite  and  arsenopyrite 
are  the  earliest  phases  followed  by  sphalerite,  chalcopyrite ,  galena, 
stibnite  and  sulfosalts  forming  in  the  late  stages. 

Ridland  (1941)  related  the  whole  process  of  ore  deposition  to  a 
single  ore-bearing  solution,  introduced  into  the  shear  zone,  from  which 


FIGURE  3 


Paragenetic  Relationships  of  the  Ore  and 
Gangue  Minerals  from  Con  Mine,  Yellowknife 
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all  the  ore  minerals  were  formed,  whereas  Graham  and  Kaiman  (1951)  out¬ 
lined  two  main  periods  of  mineral  formation,  separated  by  a  period  of 
fracturing.  The  authors  did  not  mention  whether  the  two  phases  were  related 
to  two  different  generations  of  ore  solutions  or  not.  Coleman  (1957) 
envisioned  three  distinctly  different  and  separate  stages  of  mineralization, 
each  involving  the  introduction  of  a  solution  having  a  composition  different 
from  those  which  preceded  them,  or  two  main  stages,  with  a  sub-phase 
during  the  second  stage.  He  explained  the  manner  of  emplacement  to  be  by 
filling  susceptible  zones  by  ore-fluids,  accompanying  deformation. 

During  the  present  study  a  suite  of  samples  from  the  Campbell  shear 
zone  were  examined.  Particular  attention  was  paid  to  the  textural 
relationships  between  the  ubiquitous  sulfides,  pyrite,  galena,  sphalerite 
and  chalcopyrite .  The  results  of  these  investigations  have  supported  the 
previous  work  of  Ridland  (1941) ,  Graham  and  Kaiman  (1951)  and  Kaiman 
(1951).  The  paragenetic  relationship  of  the  sulfide  minerals,  gold  and 
gangue  minerals  is  diagramatically  represented  in  Fig.  3  and  pertinent 
illustrations  and  descriptions  follow. 

Pyrite  and  arsenopyrite ,  along  with  galena  and  sphalerite  apparently 
belong  to  the  primary  phase  of  mineralization.  Pyrite,  like  quartz  is 
ubiquitous  in  the  system  and  appears  to  have  continued  to  be  deposited 
during  all  phases  of  ore-deposition.  Galena  and  sphalerite  were  probably 
deposited  during  two  distinct  phases  of  mineralization.  The  very  character¬ 
istic  feature  noticed  in  some  of  the  samples  is  the  presence  of  minute 
crystals  of  pyrite  enclosed  and  oriented  parallel  to  the  vein  pattern  in 
which  sphalerite  is  the  major  phase  (matrix  mineral).  Chains  of  pyrite 
crystals  parallel  the  boundaries,  as  well  as  the  centre  parts  of  the 
sphalerite  veinlets  in  quartz  (Fig.  4a).  Various  sulfides  present  in  the 


FIGURE  4  (a) 


Euhedral  crystals  of  pyrite  isolated  and  parallel  to  the 
walls  of  the  sphalerite  vein  (light  grey) .  Quartz  (dark 
grey)  is  present  on  both  sides  of  the  vein 


FIGURE  4  (b) 

Replacement  vein  of  galena  (light  grey)  and  chalcopyrite 
(dark  grey)  in  pyrite.  Galena  and  pyrite  show  mutual  mineral 


boundary  relations 


I 


I 


FIGURE  5 


Pyrite  and  sphalerite  -  rhythmic  banding  associated  with 

quartz  (dark  grey) 
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hydrothermal  solutions  might  attain  the  necessary  degree  of  saturation 
thus  resulting  in  the  precipitation  of  that  particular  phase.  It  is 
thought  that  even  though  all  the  sulfides  may  be  contemporaneous,  it  is 
probable  that  the  parochial  precipitation  of  sulfides  in  certain  sectors 
might  have  been  influenced  by  local  physico-chemical  factors.  As  a 
result,  replacement  textures  and  vein  deposition  of  other  sulfides  in 
pyrite  grains  might  occur  (Fig.  4b).  Some  conflicting  mutual  relation¬ 
ships  were  also  noted.  However,  in  the  opinion  of  Coleman  (1957)  the 
minerals  formed  at  any  single  period  of  mineralization  exhibit,  for  the 
most  part,  conflicting  mutual  relationships  or  else  relationships  that 
suggest  contemporaneity.  The  terms  ’time  of  mineralization',  ’early’  and 
’late’  are  strictly  used  in  a  relative  sense. 

Successive  deposition  inferred  by  crustif ication  of  two  or  more 
minerals  has  been  interpreted  by  Edwards  (1965)  as  being  due  to  fractional 
crystallization  and  could  be  regarded  as  a  special  varient  of  simultaneous 
crystallization.  In  Fig.  5,  it  is  shown  that  pyrite  and  sphalerite  show 
quite  distinct  alternate  banding  with  clear-cut  boundaries.  This  is 
presumably  because  of  different  salts  in  aqueous  solution  attaining 
saturation  with  a  rhythmic  periodicity.  In  Bastin' s  opinion  (1931),  it 
is  not  possible  to  envisage  solutions  changing  their  compositions  radi¬ 
cally  for  the  deposition  of  each  layer,  but  it  could  be  reasonably  inter¬ 
preted  as  being  caused  by  delicately  changing  conditions  of  individual 
saturations.  Thus  crustif ication ,  probably,  represents  local  modifica¬ 
tions  of  physico-chemical  factors  as  mentioned  earlier. 

Sphalerite  was  observed  showing  local  variations  in  colour  at  differ¬ 
ent  elevations.  The  type  of  sphalerite  sampled  at  3512  E  displays  a  dark 
reddish  to  dark  yellow  colour,  indicating  higher  percentages  of  Fe,  where 
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as  light  yellow  to  honey  yellow  low  Fe  sphalerite  was  collected  at  the 
3300  ft  level.  Further  work  obviously  needs  to  be  performed  on  the 
study  of  this  mineral  from  the  mine  to  determine  the  physico-chemical 
controls  on  its  composition. 

3.3  PARAGENETIC  RELATIONSHIP  OF  GOLD  TO  OTHER  MINERALS 

The  association  of  gold  with  vein-quartz  is  so  common  and  well 
known  that  no  detailed  textural  studies  were  attempted.  The  gold  ores 
of  the  Yellowknife  area  have  been  extensively  studied  by  various  authors 
(e.g.,  Coleman  1953)  and  diverse  opinions  have  been  forwarded  regarding 
the  paragenetic  position  of  gold  deposition.  Coleman  (1953)  stated 
that  the  gold  is  mainly  associated  with  sulfosalts  of  a  late  stage  in 
the  paragenetic  sequence.  Ridland  (1941)  stated  that  the  gold  is  commonly 
associated  with  the  metallic  minerals  and  quartz  which  were  introduced 
during  the  first  stages  of  ore  deposition,  whereas  the  carbonates  were 
introduced  during  the  second  stage.  A  long  period  of  gold  deposition 
is  postulated  by  Ridland  (1941) ,  as  the  gold  occurs  in  two  different 
forms,  namely  (a)  as  finely  segregated  particles  distributed  within 
fractures  or  along  grain  boundaries  and  (b)  associated  with  sulfosalts. 
Another  mode  of  occurrence  of  the  gold  was  reported  by  Coleman  (1953) 
who  noted  the  gold  occurring  as  minute  blebs  in  homogenous  pyrite, 
sphalerite  or  arsenopyrite .  These  blebs  do  not  exhibit  any  connection 
with  the  minute  fractures  in  the  quartz.  Probably  this  might  be  the 
gold  cogenetic  with  arsenopyrite.  However  Coleman  (1953)  believes 
that  the  earlier  coarse-grained  gold  might  have  been  partly  remobilized 
and  redeposited  by  later  mineralizing  fluids  as  auro  complexes. 

It  is  noteworthy  that  at  Kolar,  India,  much  of  the  gold  occurs  within 


narrow  veins  associated  with  other  minerals  in  quartz,  which  under¬ 
went  fracturing  and  recrystallization  prior  to  the  advent  of  gold 
(Pryor ,  1924) . 
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CHAPTER  IV 

EXPERIMENTAL  WORK 

4.3.  SULPHUR  ISOTOPE  STUDIES: 

(a)  Introduction 

Studies  of  variations  in  sulfur  isotope  distribution  within  an 
ore-body  may  be  employed  to  elucidate  the  physico-chemical  factors 
influencing  ore  deposits  and  consequently  provide  clues  as  to  the 
genetic  relationships  of  the  ores  (Sakai  1968) .  During  ore  transport 
and  deposition,  isotope  exchange  can  be  visualized  to  have  been  oper¬ 
ative  between  coprecipitating  sulfide  minerals,  as  well  as  between  the 
precipitates  and  the  sulfide  ions  in  the  surrounding  fluids.  Recently 
sulfur  isotope  fractionation  theories  have  been  utilized  in  geothermo¬ 
metric  studies  (Kajiwara  et  al.,  1970;  Rye,  1970)  and  in  the  investi¬ 
gation  of  the  conditions  of  deposition  of  various  sulfides  and  the 
nature  and  isotopic  composition  of  the  hydrothermal  fluids. 

(b)  Experimental  procedure: 

The  sulfide  phases  for  analysis  were  separated  by  crushing,  taking 
the  170-230  mesh  fraction,  and  separating  the  various  components  utilizing 
a  Franz  isodynamic  separator,  heavy  liquid  methods  (tetrabromethane)  and 
final  purification  by  hand  picking. 

A  10  mg  sample  of  the  sulfide  was  then  mixed  thoroughly  with  100  mg 
of  CuO  (pre-treated)  and  inserted  into  a  quartz  tube  with  quartz  wool 
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loosely  packed  on  both  sides  of  the  sample.  This  sample  was  then 
heated  in  vacuum  at  1040°C  for  at  least  10  minutes  and  the  evolved 
H^O,  02  and  C02  fractions  were  then  removed  by  using  freezing  traps  of 
liquid  air,  and  a  mixture  of  acetone  and  dry-ice.  The  S02  analysis  line 
is  illustrated  in  Fig.  6.  The  pure  So2  was  then  trapped  into  a  standard 
break-seal. 

The  sulfur  isotope  analyses  were  performed  on  a  12",  90”  radius 

34 

Neir-type  Mass-spectrometer  employing  the  YKS-GMIA  pyrrhotite  of  S  -1.94%0 
as  the  laboratory  standard  (relative  to  Canon  Diablo),  by  introducing 
the  standard  and  the  sample  alternately  and  collecting  masses  64  and  66. 

The  ratios  of  these  masses  were  obtained  by  direct  printout. 

The  results  of  the  isotopic  analyses  are  represented  in  the  usual 
permil  notation  defined  by: 


s34%.  = 


(s34/s32)  -  (S34/S32) 

/b  •’sample  1  ;  standard 


(s34/s32) 


standard 


x  1000 


by  applying  the  following  corrections:  (a)  a  correction  factor  ?B  against 

the  working  standard  which  was  derived  to  be  0.0 5  in  this  case,  (b)  a 

fractionation  correction  for  the  standard  gas  with  time  (0.99820),  and 

(c)  a  correction  factor  of  oxygen  which  was  taken  to  be  1.0907. 

Certain  samples  were  tested  for  reproducibility  by  replicate  analysis 

and  the  range  obtained  was  ^-0.07  to  0.18%o.  For  the  purposes  of  compar- 

o/  22 

ison  with  other  published  S  /S  values,  all  the  values  obtained  were 

recalculated  with  respect  to  the  common  standard  of  troilite  from  the 

32  34 

Canon  Diablo  meteorite,  which  has  a  ratio  of  S  /S  =  22.21%0or 
(S^/S32  =  0.0/Q.  The  standard  derivation  of  the  results  was  calculated 
individually  for  all  the  samples.  For  90%  of  the  samples  the  standard 


FIGURE  6 


(in 


Schematic  Diagram  representing 
Sc>2  Extraction  Train 
the  Geology  Department,  University 


the 

of  Alberta) 


CM 


O 

C/0 


P, 

§ 

P. 

P 

a 

CO 

> 


■r) 


cn 

pu 

cO 

u 

H 


29 


deviation  falls  in  the  range  of  0.01  to  0.07%o  whereas  for  the  rest,  the 
standard  deviation  ranges  from  0.09  to  0.12. 

(c)  Results 

The  analytical  results  of  the  isotopic  studies  on  various  sulfide 
minerals  from  the  Con  Mine  together  with  the  descriptions  of  those  minerals 
are  presented  in  Table  2.  The  sample  locations  also  being  indicated  on  the 
longitudinal  section  Fig.  7.  The  results  are  diagrammatically  represented 
in  Fig.  8.  In  the  table,  depth  indicates  the  stopes  or  drift  at  a  particular 
elevation  represented  by  the  figure. 

4.2  FLUID  INCLUSION  STUDIES 

(a)  Introduction 

Homogenization  temperature  determinations  on  the  fluid  inclusions  in 
quartz  and  calcite  have  long  been  utilized  to  confirm  the  data  obtained  from 
other  methods  of  geothermometry  and  to  discuss  some  of  the  related  physical 
parameters  in  the  light  of  the  fluid  phase  and  gas  phases  present.  The 
validity  of  the  data  obtained  from  fluid  inclusion  studies  has  long  been 
debated  on  the  basis  of  the  fact  that  the  filling  temperatures  obtained  may 
not  actually  represent  the  temperatures  of  ore  deposition,  even  after  apply¬ 
ing  the  pressure  corrections.  In  addition  to  that,  the  analyses  of  the  ex¬ 
tracted  fluids  from  the  fluid  inclusions  were  doubted  to  be  the  represent¬ 
atives  of  the  true  composition  of  the  soluble  salts  within  the  ore  solutions. 
However,  the  results  of  considerable  experimental  work  concerning  epithermal 
deposits  in  Russia  (Borina  1963)  exhibit  a  high  degree  of  consistence,  in¬ 
dicating  that  the  above-mentioned  determinations  are  quite  meaningful. 

(b)  Homogenization  Temperatures 

The  filling  temperatures  were  determined  on  a  heating  stage  designed 
by  Ohmoto.  Few  samples  were  found  to  be  suitable  for  filling  temperature 
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TABLE  2 


Analytical  Results  and  Sample  Descriptions 
(Sulphides  mainly  Pyrite,  Galena,  Sphalerite 
and  Chalcopyrite  from  Con  Mine,  Yellowknife) 


SAMPLE 

NO 

MINERAL 

s34% 

LOCATION 

DEPTH  (LAT) 

DESCRIPTION 

1. 

Pyrite 

+1.5 

2927 

Ay-1 

(19200) 

Py  with  quartz,  chip  sample 

2. 

Pyrite 

+0.7 

3512 

E  -1 

(19700) 

Py  with  other  sulfides  +  qutz 

3. 

Pyrite 

+1.2 

4111 

Ax-1 

(19100) 

Py  +  quartz  +  sulfides 

4. 

Pyrite 

+1.4 

4111 

Ax-2 

(19100) 

Py  with  quartz 

5. 

Pyrite 

+1.3 

4311 

Ax-2 

(19050) 

Py  with  other  sulf  +  quartz 

6. 

Sphalerite 

-0.5 

3322 

Ay-1 

(19600) 

Sph,  gal  and  quartz 

7. 

Sphalerite 

-0-6 

4311 

Ax-1 

(19050) 

Sph  with  quartz 

8. 

Galena 

-1.6 

3322 

Ay-1 

(19600) 

Galena  with  quartz 

9 

Galena 

-2.6 

3512 

E  -1 

(19700) 

Galena  with  other  sulf  +  qutz 

10. 

Galena 

-2.8 

4311 

Ax-1 

(19050) 

Gal  and  other  sulf  +  qutz 

11. 

Pyrite 

+0 . 4 

4511 

DS-1 

(18900) 

Pyrite  +  quartz 

12 

Sphalerite 

+0.9 

4311 

Ax-4 

(19100) 

Pyrite  +  sph  +  quartz 

13. 

Pyrite 

+1 . 6 

3512 

A  -1 

(19450) 

Pyrite  +  quartz 

14. 

Pyrite 

+1 . 4 

4513 

DS-1 

(18700) 

Pyrite  and  quartz 

15. 

Galena 

-3.4 

4311 

Ax-2 

(19050) 

Gal  +  quartz  asso  py 

16. 

Pyrite 

+1.1 

4311 

Ax-2 

(19050) 

Py  +  other  sulf  +  quartz 

17. 

Pyrite 

+1.2 

3512 

E  -3 

(19700) 

Py  +  other  sulf  +  quartz 

18. 

Pyrite 

+0.4 

3322 

Ay-1 

(19600) 

Pyrite  +  quartz 

19. 

Chalcopyrite 

+0.5 

3512 

E  -4 

(19700) 

Cpy  +  pyrite  +  quartz 

20. 

Pyrite 

+0 . 6 

3512 

E  -5 

(19700) 

Pyrite  +  sulf  +  quartz 

21. 

Pyrite 

+1.4 

3512 

E 

(19700) 

Py  +  massive  sulf  +  qutz 

22. 

Pyrite 

+1 . 6 

3711 

Ax-2 

(19100) 

Py  assoc  qutz  and  sulf 

23. 

Galena 

-0.7 

3711 

Ax-2 

(19100) 

Gn  assoc  py  and  quartz 

(continued. . . ) 


(continued) 


SAMPLE 

NO 

MINERAL 

s34% 

LOCATION 

DEPTH  (LAT) 

DESCRIPTION 

24. 

Pyrite 

+1.2 

4311  Ax-1 

(19050) 

Pyrite  +  quartz  and  sulf 

25. 

Chalcopyrite 

+0.4 

4311  Ax-1 

(19050) 

Cpy  assoc  with  py  and  qutz 

26. 

Galena 

-1.9 

3321  Ay-1 

(19250) 

Galena  +  qutz  +  pyrite 

27. 

Sphalerite 

+1.0 

4311  Ax-T 

(19100) 

Sph  +  qutz  +  sulf  minor 

28. 

Sphalerite 

+0.9 

3511  H  -T 

(20140) 

Sph  +  qutz  +  minor  sulf 

29. 

Pyrite 

+1.9 

3512  A  -T 

(19440) 

Pyrite  +  quartz 

30. 

Pyrite 

+0.4 

4115  Av-T 

(18830) 

Pyrite  +  quartz 

31. 

Pyrite 

+1.6 

3711  Ax-1 

(19100) 

Pyrite  +  quartz 

32. 

Galena 

-2.7 

3322  Ay-1 

(19600) 

Galena  with  quartz 

33. 

Galena 

-2.0 

3511  H  -T 

(20140) 

Gal  +  qutz  +  minor  sph 

34. 

Galena 

-3.0 

4115  Av-T 

(18830) 

Gal  +  qutz  +  minor  py 

35. 

Pyrite 

+1.7 

4311  Ax-T 

(19100) 

Pyrite  +  quartz  +  sulf 

36. 

Pyrite 

+1 . 3 

2927  Ay-T 

(19140) 

Pyrite  +  quartz 

37. 

Galena 

-0.8 

3321  Ay-1 

(19000) 

Gal  +  other  sulf  +  qutz 

38. 

Pyrite 

+1.4 

4311  Ax-4 

(19050) 

Py  +  other  sulf  +  qutz 

39. 

Sphalerite 

+0.5 

4311  Ax-4 

(19050) 

Sph  +  qutz  +  minor  sulf 

40. 

Pyrite 

+1 . 6 

3512  A 

(19450) 

Pyrite  +  quartz 

41. 

Pyrite 

+1.3 

3312  A 

(19700) 

Pyrite  +  quartz 

42. 

Pyrite 

+1.7 

3711  Ax-2 

(19100) 

Py  +  other  sulfides 

43. 

Pyrite 

+1.7 

3511  H  T 

(20140) 

Py  +  qutz  +  minor  sulf 

FIGURE  7 


Longitudinal  Section  Campbell  Shear  System 
Con  Mine,  Yellowknife 

(The  first  two  digits  refer  to  the  depth  in  1000’ s  of  feet 
and  the  rest  the  S  tope  location  with  respect  to 
foot  wall  -  hanging  wall  relationship) 


(map  enclosed  in  the  pocket) 


FIGURE  8 


Distribution  of  Sulfur  Isotope  Values 
in  different  Sulfide  Minerals, 

Con  Mines,  Yellowknife 
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determinations  as  most  of  the  fluid  inclusions  were  almost  beyond  the 
resolution  of  the  microscope  and  therefore  gas-liquid  phases  could  not 
be  distinguished.  Furthermore,  since  the  quartz  is  so  often  strained, 
polished  sections  tend  to  be  rather  opalescent  and  translucent. 

However,  four  of  the  fluid  inclusions  observed  in  quartz  showed  filling 
temperatures  of  270-10°C  and  the  inferred  temperatures  of  the  ore 
solutions,  when  corrected  fall  in  the  range  of  300-350°C  for  an  average 
pressure  of  400-500  atmospheres  (Lemmlein  and  Kletsov,  1961) . 

4.3  INVESTIGATIONS  OF  K/Na  RATIOS  IN  FLUID  INCLUSIONS 
(a)  Introduction 

Sawkin  (1966)  indicated  that  there  is  a  relative  variation  of  K/Na 
ratios  in  natural  waters  with  temperature.  The  empirical  plot  of  White 
(1968),  when  Gammon's  experimental  data  is  superimposed,  also  showed  a 
relationship  between  temperature  and  K/Na  ratios.  An  increase  in  the 
temperature  of  chloride  water  in  contact  with  feldspar  and  micas  will 
apparently  shift  the  fluid-mineral  equilibria,  thereby  increasing  the 
K/Na  ratio  in  the  fluid  phase,  as  Na  is  exchanged  for  K  in  the  minerals 
(Orville  1963  Hemley  1967).  When  hydrothermal  fluids  evolved  from 
granite  melts,  react  with  the  country  rock,  exchange  type  reactions  can 
be  expected,  and  equilibrium  is  obtained  as  a  result  of  the  formation  of 
alteration  zones  in  the  form  of  aureoles  surrounding  the  veins.  In  general 
the  Na/K  ratios  for  solutions  in  contact  with  rhyolitic  volcanic  rocks 
during  hydrothermal  alteration  were  found  to  range  from  10  to  15  at  250°C 
to  approximately  4  at  500°C  (Ellis  and  Mahon,  1967).  Thus  it  can  be 
shown  that  there  is  an  influence  of  temperature  on  the  K/Na  ratios  in  the 
ore-fluids.  Furthermore,  K/Na  ratios  would  help  to  elucidate  the 
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depositional  history,  temperature  gradients  and  equilibrium  relations 
which  will  be  discussed  in  detail  later  in  this  thesis. 

(b)  Experimental  procedures 

K/Na  ratios  of  the  fluid  inclusions  in  quartz  from  the  Con  Mine 
were  determined  using  a  leaching  procedure  after  Sawkins  (1966)  and 
Grooves  and  Soloman  (1969) .  In  order  to  determine  the  relative  amounts 
of  K  and  Na  in  the  fluid  inclusions  from  the  samples  collected  at 
different  elevations  in  the  Campbell  shear  system,  six  quartz  samples 
were  carefully  cleaned  and  crushed  separately  in  a  clean  tungsten- 
carbide  mixer  mill  for  four  minutes.  The  crushed  material  was  leached 
with  distilled  and  demineralized  water,  washing  the  mixer  about  six 
to  eight  times  for  every  sample.  Final  washings  were  performed  with 
1.0  mole  nitric  acid.  The  resulting  slurry  was  kept  overnight  and 
filtered.  The  filtrate  was  heated  gently  to  dryness.  The  volume  of 
all  the  samples  as  well  as  the  two  blank  runs,  were  adjusted  to  25  cc 
and  K  and  Na  concentrations  of  these  solutions  were  determined  on  the 
flame-photometer . 

The  experimental  accuracy  of  these  determinations  was  determined 
by  two  blank  runs,  namely  (a)  following  the  same  procedure  with 
demineralized  water  without  any  sample,  and  (b)  by  following  the  same 
procedure  using  synthetic  quartz  as  blank.  The  results  obtained  should 
be  viewed  in  the  light  of  the  following  limitations. 

(i)  The  accuracy  of  the  results  is  difficult  to  estimate 

owing  to  the  ubiquity  of  secondary  inclusions  in  all  samples. 


The  higher  concentrations  of  K  in  these  secondary 
inclusions  might  be  expected  to  be  due  to  latter 
metamorphic  and  metasomatic  effects,  however.  Such 
conclusions  are  merely  speculative  unless  certain 
geological  evidences  are  shown.  No  specific  limitations 
are  possible  to  impose  on  the  extent  of  contamination 
by  this  factor  at  this  stage. 

(ii)  The  number  of  fluid  inclusions  in  the  quartz  samples 
varied  considerably  from  one  sample  to  the  other. 

Thus  the  presence  of  variable  concentrations  of 
secondary  inclusions  in  the  same  samples  causes  the 
results  obtained  to  show  a  lack  of  internal  consistency. 
Extreme  caution  was  taken  throughout  the  experimental 
procedure,  however,  unanticipated  contamination  might 
have  been  present  in  any  one  of  the  individual 
samples . 

(iii)  Even  though  most  of  the  quartz  was  hand-picked  for  this 
purpose,  the  quartz  might  still  (in  certain  cases) 
contain  muscovite  as  very  fine-grained  inclusions;  (some 
of  the  samples  contained  fine  muscovite  associated  with 
quartz.  Some  of  the  samples  contained  fine  muscovite 
associated  with  quartz  and  this  contamination  is  believed 
to  have  contributed  a  major  proportion  ion  the  higher  K/Na 


ratios  obtained) . 
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(c)  Results 

The  results  of  the  determinations  are  tabulated  in  Table  3. 

To  obtain  an  overall  value  of  K/Na  ratio,  a  plot  is  drawn  between  K 
and  Na  (ppm)  (Fig.  9)  and  it  appears  that  some  correlation  exists 
between  the  K  concentration  and  the  corresponding  Na  concentration 
in  any  given  sample.  All  the  values  have  been  corrected  for  blanks. 

The  final  values  obta-ned  from  the  graph  is  0.66^0. 05  (atomic 
ratio) . 

4.4  OXYGEN-CARBON  ISOTOPE  STUDIES  -  CARBONATES 

(a)  Introduction 

Oxygen  isotope  fractionation  studies  have  been  applied  to  the 

study  of  hydrothermal  processes  and  facilitated  the  determination  of 

temperatures  of  formation  of  minerals  and  the  determination  of  the 

origin  and  source  of  the  solutions  involved.  Furthermore,  in  the  light 

of  oxygen  isotope  fractionation  studies,  the  validity  of  the  assumption 

of  equilibrium  among  the  mineral  species,  has  been  discussed.  It  has 

18  16 

also  been  demonstrated  that  the  0/0  of  water  in  equilibrium  with 
the  mineral  assemblages  at  a  given  temperature  could  be  calculated 
by  applying  the  experimentally  determined  mineral-water  relations 
(Clayton  1969)  . 


c 


TABLE  3 


Concentration  of  K  and  Na  in  the  Leached  Solutions 
from  Quartz  Samples,  Con  Mine,  Yellowknife 
(all  the  readings  are  corrected  for  the  blank) 


NO 

SAMPLE 

SPECIFICATION 

K  (ppm) 

Na  (ppm) 

K/Na 

(at  ratio) 

Na/K 
(at  rat 

1. 

2927  Ay 

27.0 

26.5 

0.599 

1.67 

2. 

3322  Ay 

45.0 

31.0 

0.853 

1.17 

3. 

3512  E 

33.0 

39.5 

0.491 

2.03 

4. 

3711  Ax 

13.5 

11.0 

0.721 

1.39 

5. 

4111  Ax 

41.0 

40.0 

0.602 

1 . 66 

6. 

4511  DS 

26.5 

23.5 

0.663 

1.51 

7. 

Blank  1 

0.0 

0.0 

8. 

Blank  2 

3.5 

4.5 
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FIGURE  9 


Plot  showing  the  Relationship  between  the  K  and  Na  (experimental  data) 
determined  by  flame  photometry,  in  the  leached  solutions  from  quartz 
samples  from  the  Con  Mine.  The  general  trend  of  K  and  Na  concentration 
is  shown  by  the  straight  line. 

The  K+/Na+  (atomic  ratio)  from  the  plot  is  0.66 


DISTRIBUTION  OF  K  -  Na  (in  ppm.)  IN  FLUID  INCLUSIONS  OF  QUARTZ 


t 


Note : 


to  obtain  K/Na  atomic  ratio, 
weight  ratio  by  0.588. 


multiply 
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(b)  Results 

Six  carbonate  samples  were  collected  from  various  parts  of  the 
mine  workings  and  Co^  was  extracted  by  using  an  oxygen  line  for 
carbonates  designed  and  constructed  by  P.  Fritz  at  the  University  of 
Alberta.  The  carbondioxide  gas  was  analyzed  on  a  1/2",  90°  radius 
Nier-type,  mass  spectrometer  using  a  laboratory  standard  of  Fisher 
scientific  company  calcium  carbonate.  The  following  results  for  the 
isotopic  composition  of  oxygen  and  carbon  in  the  hydrothermal  carbonates 
were  obtained  after  correcting  for  (a)  tail  effects,  (b)  background,  and 

(c)  mass  contribution  effects.  All  the  values  are  reported  in  the 
standard  %.  permil  notation  with  reference  to  the  standards  mentioned 
below: 


(Against  SMOW  for  Oxygen  and  PDB  for  Carbon) 


SNO 

Prep.  No. 

Collection  (locality) 

018%o 

C13%, 

1 

1453 

3511  St  ReX^  calcite 

+10.6 

-4.2 

2 

1284 

2900  Drift  calcite 

+8.4 

-3.6 

3 

1376 

3511  St  X^  calcite 

+9.6 

-0.5 

4 

2900 

2900  Drift  X^  calcite 

+8.7 

-2.9 

The  results  are  reported  in  the  familiar  notation,  where 

R  R 

.  .18,,.  ,  sample  -  standard, 

SO  %.  =  [  - 1 - - - ]  X  1000 

^standard 
18  16  13  12 

where  R  represents  the  0/0  or  C  /C  ratio  of  Co^  from  carbonate 
samples.  The  standards  recalculated  for  are  SMOW  for  oxygen  isotopes 


and  PDB  for  carbon  isotopes. 
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CHAPTER  V 

INTERPRETATIONS  -  DISCUSSION 

5.1  GEOTHERMOMETRY  -  SULFUR  ISOTOPES 

Isotopic  fractionation  of  sulfur  between  synthetic  sulfide 
mineral  pairs  (FeS^-PbS,  ZnS-PbS  and  FeS^-Zns)  as  a  function  of 
temperature  and  duration  of  reaction  time  has  been  investigated  by 
many  authors  (Kajiwara  1969,  1970;  Rye  1969;  Grootenboer  1969).  From 
the  results  shown  in  the  table  2,  the  isotopic  composition  of  coexisting 
sulfide  pairs  have  been  utilized  to  calculate  the  per  mil  deviations 
between  the  specific  sulfide  pairs;  the  computations  being  given  in  the 
table  4.  With  respect  to  the  selection  of  the  appropriate  geothermo¬ 
metric  plots,  shown  in  Fig.  10,  all  these  plots  were  taken  into  consider¬ 
ation  in  the  initial  stages  and  subsequently  comparative  values  were 
obtained  by  employing  other  geothermometric  procedures,  such  as  fluid- 
inclusion  filling  temperature  methods.  Experimental  and  theoretical 
plots  are  presented  in  Fig.  10  and  11. 

The  nature  of  the  distribution  of  isotopic  composition  of  sulfur  in 
pyrite,  sphalerite,  chalcopyrite  and  galena  shows  a  general  fractionation 
trend  as  theoretically  calculated  by  Sakai  (1968)  based  upon  the  partition 
ratios  for  crystalline  FeS^,  PbS  and  ZnS.  The  theoretical  trend  of 

fractionation  of  sulfur  isotopes  in  coexisting,  cogenetic  sulfide  species, 

34 

as  postulated  by  Sakai  (1968),  is  that  the  S  enrichment  would  be: 
pyrite  >  pyrrhotite  _>  sphalerite  >  chalcopyrite  )>  galena.  This  same 


TABLE  4 


Temperature  determinations  based  upon 
of  sulfur-isotope  fractionation  among 
Rye  1969,  and  theoretical  plots  Sakai 


the  experimental  plots 
sulfides  (Kajiwara  1970; 
1968;  Houlston  1958) 


Kg  =  (S34/S32)  /  (S3Vs32) 
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N.B.  1.  Sample  locations  refer  to  the  specified  working  stopes  on  the  levels,  represented  by  the  left 
two  digits  in  thousands  of  feet  depth. 

2.  Numbers  specified  under  ’sample  pair'  represent  the  same  in  Table  3  (isotopic  compositions  of 
sulfides) . 

3.  Py-pyrite;  Sph-sphalerite ;  Gn-galena. 


FIGURE  11 


Theoretical  sulfur-isotope  fractionation  curves  for  pyrite- 
galena,  pyrite-sphalerite ,  sphalerite-galena ,  based  upon 
Sakai's  (1968)  estimates  of  partition  function  ratios  for 
crystalline  FeS^ ,  PbS,  ZnS.  is  approximately  equal  to 
1000  In  ^ and  for  the  exchange  reactions  given,  is  identical 
to  the  equilibrium  constant  K. 

Isotope  data  obtained  from  co-existing  sulfides  (Con  Mine, 
Yellowknife)  are  also  presented  on  the  plot. 
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trend  is  clearly  observed  in  the  samples  collected  from  the  Campbell 
shear  system.  The  trend  might  imply  that  the  mineral  pairs  have 
attained  equilibrium  or  closely  approached  it.  However,  sphalerite  and 
chalcopyrite  show  a  slight  variation  which  could  be  either  due  to  local 
isotopic  variations  or  due  perhaps  to  experimental  errors  (Fig.  12). 

A  close  observation  of  the  spread  of  sulfur  isotope  composition  of 
the  three  major  sulfides  (Fig.  9),  namely  pyrite,  galena  and  sphalerite, 
together  with  corresponding  textural  studies  on  these  three  mineral 
pairs,  indicate  a  probable  history  of  precipitation  of  the  sulfides 
which  may  be  as  follows:  pyrite  was  ubiquitous  within  the  shear  system 
and  no  temperature  gradient  was  attained  during  the  precipitation  of 
the  sulfides.  However,  other  physico-chemical  factors,  such  as  the 
hydrogen-ion  concentration  and  fugacity  of  oxygen,  might  have  varied 
considerably,  effecting  changes  in  the  isotopic  composition  of  the  sulfur 
species  within  the  aqueous  solutions.  This  has  been  treated  theoretically 
by  Ohomoto  (personal  communication)  and  is  to  be  discussed  in  detail 
later . 

The  galena  and  sphalerite  exhibit  features  of  simultaneous  precipit- 

34 

ation  along  with  pyrite  and  the  corresponding  S  values  of  galena 
and  sphalerite  show  a  range  of  -2.9%0  to  -2.0%o  and  -0.6 %c  to  -0.5%o  re¬ 
spectively,  while  pyrite  showing  an  isotopic  composition  ranging  from  +0.5 
to  1.0%o.  However,  where  there  are  veins  and  replacement  features,  namely 
sphalerite  and  chalcopyrite  replacing  pyrite,  and  galena  replacing  pyrite, 
the  isotopic  compositions  of  sulfur  of  these  sulfides  are  distinctly 
heavier,  relative  to  those  described  earlier.  The  values  of  the  isotopic 
compositions  of  sulfur  in  galena  and  of  sphalerite  would  have  been 
considerably  modified  during  the  second  phase  of  precipitation  indicated 


£  ■ 


FIGURE  12 


Plot  showing  the  mutual  relationship  between  the  pyrite- 
sphalerite  and  pyrite-galena  (theoretical  data  of  Sakai 
1968)  and  superimposed  Con  Mine  data. 
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by  the  textural  studies.  Likewise  the  pyrite  also  shows  relatively 
34 

heavier  S  values  ranging  from  1.3  to  1.8 %0.  When  the  earlier 
sulfide  pairs  are  accepted  as  having  been  precipitated  simultaneously, 
the  temperature  determinations  on  these  pairs  would  reveal  the  temper¬ 
ature  of  formation  of  cogenetic  sulfide  pairs  if  the  other  factors  are 
the  same  during  the  precipitation  of  all  the  sulfides.  The  highest 
temperatures  determined  on  the  pairs  range  from  340°-360°  whereas  the 
lowest  values  could  be  considered  to  be  within  the  range  of  310°-300°, 
a  temperature  range  which  is  apparently  quite  common  in  hydrothermal 
deposits.  The  range  of  temperatures  determined  for  the  Con  Mine  deposit 
is  between  300-350°C,  i.e.,  the  temperature  interval  might  indicate  the 
interval  during  which  the  sulfides  and  other  minerals  have  been  precipi¬ 
tated  from  the  hydrothermal  solutions.  The  temperatures  obtained  for  the 
ore  deposits  of  the  Campbell  shear  zone  seem  to  be  quite  meaningful  in 
the  light  of  the  earlier  studies  (Kaiman  et  al.  1964). 

As  the  precipitation  of  pyrite  continues,  toward  the  latter  phase  of 
the  depositional  episode,  sphalerite  and  galena  would  have  reached  super¬ 
saturation  due  to  the  changing  physico-chemical  characteristics,  and 
hence  had  developed  heavier  characteristics,  compared  with  the  previous 
galena  and  sphalerite.  The  isotopic  composition  of  galena  and  sphalerite 
range  from  -1.5%0  to  -1.0%o  and  0.4%o  to  0.9%o  respectively,  while  pyrite 
shows  an  isotopic  composition  ranging  from  1.0%o  to  1.7%0.  This  is 
illustrated  schematically  in  Fig.  13. 

The  geochemistry  and  the  origin  of  the  Yellowknife  belt  has  been 
investigated  and  critically  discussed  by  Boyle  (1961) ,  who  interpreted 
the  quartz  and  gold  as  being  epigenetic  in  character  and  being  formed 
from  silica  released  during  mineralogical  and  chemical  changes  occurring 
within  the  country  rocks  during  the  development  of  the  schist  zones. 


FIGURE  13 


Schematic  diagram  showing  the  relationship  between  the  sulfur 
isotope  variations  and  paragenetic  sequence  in  pyrite, 
sphalerite  and  galena  (Con  Mine,  Yellowknife) . 
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The  sulfides  in  the  shear  zones  are  thought  by  Boyle  (1961)  to  have 
originated  by  redistribution  of  sulfo-phile  elements  during  metamorphism 
of  the  basic  greenstone  suite  and  partly  as  the  result  of  chemical 
changes  which  took  place  during  the  development  of  schist  zones.  The 
elements  present  in  the  shear  zone  were  therefore  interpreted  as  being 
derived  from  the  country  rocks  and  were  supposedly  concentrated  by  the 
interaction  of  metamorphic  processes  within  dilatent  structures  in  the 
shear  zones. 

Wanless,  Boyle  and  Lowden  (1960)  studied  the  sulfur  distribution 
in  different  phases  of  the  lithological  units.  Their  work  is  summarized 
in  the  table  5  and  diagrammatically  represented  elsewhere  (Fig.  14  and 
15).  The  following  conclusions  were  drawn  from  their  results:  the 
available  sulfur-isotope  proportions  in  the  source  area  were  assumed  to 
be  in  the  same  proportions  as  those  found  within  meteoric  troilite  and 
thus  variations  in  the  isotopic  compositions  were  interpreted  to  be  due 
to  temperature,  pressure  and  compositional  gradients. 

(a)  The  greatest  enrichment  of  the  heavier  sulfur  isotope  took 
place  where  high  temperatures  prevailed  and  the  lighter  sulfur  isotope 
migrated  to  the  relatively  cooler  portions  of  the  lithological  units. 

The  gradients  show  higher  values  when  the  greenstones  show  lighter  sulfur 
isotope  concentrations. 

(b)  Pressure  gradients  are  responsible  for  modification  of  the 

sulfur  isotope  distribution.  Sulfides  in  the  isolated  lenses  filling 

tension  fractures  in  the  epidote-amphibolite  facies  of  the  rocks  always 
34 

contain  more  S  than  in  the  sulfides  in  the  enclosing  rocks. 

(c)  When  the  quartz-gold  bodies  formed,  the  shear-zone  was  subjected 

32 

to  intense  chemical  activity  and  S  diffused  out  of  the  shear  zone  into 
the  associated  country  rock. 
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TABLE  5 


Sulfur  isotopic  compositions  of  pyrite  collected  and  analyzed 
by  Wanless  et  al.  (1960)  from  various  lithological  units 
(adopted  from  Wanless  et  al.  1960). 


ROCK  UNIT 

NUMBER  OF 
SAMPLES 

SPREAD  OF 

S  s34%0 

AVERAGE  OF  VALUES 
(as  given  in 
Wanless  et  al.1960) 

Granitic  Rocks 

5 

+2.9  -  +6.8 

+4.0 

Greenstone  Belt 

8 

+0 . 6  —  +2 . 2 

+1.4 

Graphic  Tuffs 

6 

+2.0  -  +5.6 

+3.5 

Siliceous  Tuffs 

4 

-2.7  -  -0.3 

-1.2 

Early  Pre-quartz  lenses 
in  the  tension  fractures 

3 

+0.3-  1.4 

+0 . 8 

Con  System  Ore-deposits 

13 

+2 . 4  -  +6 . 2 

+3.8 

Negus-Rycon  System 

29 

+0 . 3  -  +3 . 6 

+1.7 

Giant  Campbell  System 

12 

+0 . 0  -  +3 . 6 

+2.0 

Parallel  Shear  Zone 

10 

0.0  -  +3 . 6 

+2.8 

FIGURE  14 


Sulfur  isotope  distribution  in  granodiorite  and  metamorphic 
facies  of  the  Yellowknife  greenstone  belt. 


FIGURE  15 

Sulfur  isotope  distribution  in  gold-quartz  lenses,  surrounding 
carbonate-sericite  schist,  chlorite  schist  and  rocks  of  epidot 


amphibolite  facies,  Yellowknife. 
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Recent  observations  by  McConnell  (1964) ,  Green  (1968) ,  and  Robertson 
and  Cumming  (1968)  and  the  author’s  discussions  and  experiences  with  the 
mine  geologists  strongly  suggest  a  hydrothermal  theory.  For  more  discuss¬ 
ions  on  this  topic  the  reader  is  directed  to  the  pertinent  literature  cited. 

5.2  PRESSURE  ESTIMATES 

It  is  extremely  difficult  to  make  a  reasonable  estimate  of  the  pressure 
of  formation  for  the  Campbell  shear  system  from  the  sparse  data  concerning 
overburden  at  the  time  of  mineralization.  The  total  pressure  in  the  system 
at  any  given  temperature  might  be  derived  by  a  number  of  independent  methods, 
namely:  (a)  from  the  plot  of  the  relationship  between  temperature  and  pressure 
at  a  given  weight  %  salt  content  in  aqueous  solutions,  obtained  by  analyzing 
the  fluid  phases  of  fluid  inclusions  for  the  Ho0-NaCl  system  (Ohmoto  and  Rye 
1970) ,  or  (b)  from  the  geobarometric  determinations  based  upon  sphalerite 
compositions  coexisting  with  pyrrhotite  or  pyrite  (Smith  1963)  or  (c)  from 
the  experimental  plots  based  upon  the  degree  of  filling  in  fluid  inclusions 
vs.  the  temperature  at  a  given  pressure  (Lemmlein  and  Klevtsov  1961). 

When  the  temperature-pressure-density  relationship  plot  of  10  wt  % 

NaCl  solution  (Ohmoto  1970)  is  considered,  the  minimum  pressure  at  350°C 
could  be  circa  150  atmospheres  which  would  represent  the  minimum  pressure 
in  the  mineralizing  system.  Boyle  (1954)  determined  the  decrepitation 
temperatures  of  fluid  inclusions  in  quartz  collected  from  the  Campbell 
shear  system  and  obtained  primary  decrepitation  temperatures  in  the 
order  of  300°C  and  secondary  decrepitation  temperatures  around  162°C. 

By  utilizing  the  filling  temperatures  derived  from  them,  in  conjunction 
with  the  plot  of  Ohmoto  (Personal  communication) ,  the  pressure 
range  obtained  is  about  850-900  atmospheres.  This  pressure  could  thus 
be  treated  as  the  maximum  at  the  time  of  deposition.  The  minimum  and 
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maximum  pressures  obtained  are  not  unrealistic,  comparing  them  with 
previous  pressure  estimates  published  for  hydrothermal  systems.  It 
is  noteworthy  that  the  application  of  pressure  corrections  and  calculation 
of  the  fugacity  of  water  vapour  are  not  greatly  affected  by  this  range  of 
pressure. 

5. 3  PRESSURE  CORRECTIONS  FOR  FILLING  TEMPERATURES 

In  fluid  inclusion  thermometry,  corrections  for  pressure  are 
absolutely  essential  in  the  calculation  of  the  correct  temperature  of 
formation  for  a  given  ore  deposit.  Diagramatic  correction-plots  have  been 
prepared  based  upon  the  experimental  work  of  previous  authors  (Lemmlein 
and  Kletsov  1961) ,  thus  enabling  the  determination  of  the  temperature  of 
formation  of  a  given  crystalline  substance.  The  temperature  corrections 
thus  obtained  for  the  Campbell  shear  zone  ranges  from  20°  to  80°C  for 
pressures  of  150  atm  and  900  atm  respectively.  When  this  correction  is 
made  for  an  average  pressure  of  500  atmospheres,  the  corrected  temperature 
would  thus  be  around  320-330°C,  whereas  the  corrected  temperatures  for 
the  maximum  pressure  of  900  atm  would  be  360°±10°C.  It  should  be  noted 
that  the  error  resulting  from  this  type  of  initial  assumption,  might  be 
within  a  range  of  1-10%.  However,  the  maximum  and  minimum  pressures  have 
been  taken  into  consideration  when  computing  the  molarities  of  the  species 
in  aqueous  solutions  and  the  effect  of  temperature  and  pressure  on 
different  sulfur  species  are  discussed. 

5.4  DEPTH  OF  ORE  DEPOSITION 

For  an  average  pressure  of  500  atmospheres,  two  extreme  limits  have 
been  considered  for  the  following  calculations.  It  was  found  that  if 
(a)  an  overhead  fluid  column  is  taken  into  consideration  (P^  =  P^)  having 
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a  density  of  0.93  gm/cc  (10  WtZ  solution  at  350°C),  and  (b)  an  average 
crustal  rock  overhead  is  considered,  the  depths  of  deposition  would  be 
in  the  order  of  5. 0-5. 5  km  and  1. 9-2.0  km  respectively  (using  the 
equation  P^_  =  pgh,  where  P^_  is  the  total  pressure  in  bars,  p  =  the  density 
of  the  column,  g  =  the  acceleration  due  to  gravity  and  h  =  the  height  of 
the  column) .  When  the  feasibility  of  these  two  values  in  the  light  of 
the  existing  geological  observations  is  considered,  the  following  points 
should  be  mentioned: 

If  the  depth  of  ore-deposition  was  2.0  km,  the  average  rock  temper¬ 
ature  at  that  depth  would  be  roughly  80°-90°C  and  the  temperatures  of 
the  ore  solutions  in  the  range  of  340°C-350°C.  Under  such  conditions 
there  would  therefore  be  a  tremendous  temperature  gradient  between  the 
ore  solutions  and  the  country  rock,  a  phenomenon  which  would  probably  be 
expressed  in  textural  and  mineralogical  features,  such  as  chilled 
marginal  zones  etc.  However,  such  features  have  never  been  observed. 

In  such  a  case  a  partial  fluid  phase  could  be  expected  in  the  hydrothermal 
plumbing  system.  An  estimate  of  the  depth  of  ore-deposition  of  4  km 
would  be  quite  appropriate  to  set  the  model  of  the  hydrothermal  system. 
Furthermore,  it  is  quite  unlikely  that  any  considerable  pressure  gradients 
could  be  expected  across  the  ore  body  except  for  very  local  pressure 
decreases  in  the  Shear  Zone. 


5.5  FUGACITY  OF  WATER  VAPOR 

The  fucacity  of  water  vapor  within  the  system  relative  to  the  total 
pressure  was  obtained  by  multiplying  the  fugacity  coefficient  by  the 
total  pressure  (Holser  1954) .  The  average  fugacity  of  water  vapor  during 
the  time  of  deposition  was  thus  determined  as  150-200  atm,  a  result  which 


is  geologically  feasible.  The  influence  of  f 


h2o 


is  discussed  in  this 
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thesis  later. 

5.6  THE  APPLICATION  OF  THERMODYNAMIC  DATA  TO  SOLUTION  CHEMISTRY 

The  application  of  thermodynamic  principles  to  solution  geochemistry 
permits  the  study  of  ore  solution  processes  and  the  evaluation  of 
various  reactions  between  electrolytes  within  the  aqueous  solutions. 

The  theoretical  approach  of  the  study  of  natural  systems  of  coexisting 
minerals  was  supported  by  Goldschmidt  (1954) .  This  author  pointed  out 
that  there  exists  a  relationship  between  the  magnitude  of  the  free 
energy  of  the  oxides,  sulfides,  sulfates  and  carbonates  of  the  various 
elements  and  their  mode  of  occurrence  in  nature.  Kordes  (1935)  indicated 
that  the  mineral  assemblages  within  the  systems  at  a  given  temperature 
could  be  mostly  influenced  by  the  fugacity  of  sulfur  in  the  systems. 

Recently  more  quantitative  studies  have  been  performed  by  Holland 
(1965)  who  has  defined  the  conditions  of  deposition  of  minerals  in 
terms  of  free  energy  and  given  physical  parameters  such  as  temperature 
and  pressure.  Holland  (1959)  applied  such  thermochemical  data  to 
oxides,  sulfides,  sulfates  and  carbonates  which  constitute  the  principal 
ore  and  gangue  minerals  of  ore  deposits.  For  any  chemical  reaction, 
the  change  in  free  energy  may  be  calculated  from  the  entropy  and  heat 
of  formation  of  the  reactants  and  products  involved  in  the  said  reaction 
and  the  equilibrium  constant  is  related  to  the  fugacity  of  ideal  gases 
in  the  system.  Thus  the  boundary  between  the  mineral  assemblages  may  be 
defined  by  equilibrium  constants  which  in  turn  are  related  to  the 
fugacity  of  sulfur  and  oxygen  in  the  system.  From  this  type  of  approach 
it  is  possible  to  specify  the  range  of  fugacities  of  the  major  gases  from 
a  knowledge  of  the  mineral  assemblages  in  any  given  hydrothermal  deposit. 
Fugacities  of  hydrogen,  carbon  monoxide,  hydrogen  sulfide  (gas)  and  SO^ 
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FIGURE  16  and  FIGURE  17 

Diagrams  showing  the  relationship  between  the  fugacity  of 
sulfur  and  oxygen  with  respect  to  t  e  mineral  stability  fields 
(specified) 

(log  f(02)  -  log  f(S2)  -  log  f(C02)  in  Fe-O-S-C,  Pb-0-  , 

Zn-O-S  systems) 


(1)  at  350°C 


figure  16  and  (2)  at  300°C 


figure  17 
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may  be  calculated  relative  to  the  known  fugacities  of  oxygen  and 

sulfur. 

5.7  DETERMINATION  OF  FUGACITIES  OF  SULFUR  AND  OXYGEN 

From  the  calculated  temperatures  of  formations,  the  fugacities  of 
oxygen  vs  fugacity  of  sulfur  diagrams  have  been  constructed  based  upon 
the  thermochemical  data  of  Holland  (1965,  1959).  From  the  major  minerals 
identified  in  the  Campbell  Shear  Zone,  the  maximum  and  minimum  fugacities 
of  oxygen  and  sulfur  have  been  fixed  as  follows  (Fig.  16  and  17) . 


Temperature 

log 

f(o2) 

log 

£(S2) 

Min 

Max 

Min 

Max 

350°C 

-31.0 

-29.0 

-10.5 

-8.5 

300°C 

-34.0 

-32.0 

-11.5 

-9.5 

The  thermochemical  data  have  been  applied  in  the  following  two  equations: 

AG  =  AT  +  BT  log  T  +  CT 
AG  «=  -RT  log  k  x  2.303 

log  k  =  AT  +  BT  log  T  +  CT  (1) 

-RT  x  2.303 

The  A,  B  and  C  values  are  reproduced  from  Holland  (1959,  1965)  in 
Appendix  1. 

The  maximum  and  minimum  pressures  have  been  considered  to  be  900 
and  150  atmospheres  to  examine  the  effect  of  pressure  upon  other 
parameters.  The  lowermost  levels  in  the  mine  reach  a  depth  of  4500  ft 
from  the  surface.  The  corresponding  corrections  for  temperature 
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determinations  on  primary  fluid  inclusions  from  quartz  samples  v/ould 
give  a  correction  range  of  1  -  10%  (which  has  been  taken  into  consideration 
in  estimating  the  probable  temperature  of  formation  of  the  fluid 
inclusions)  . 


5.8  FUGACITY  OF  CARBON  DIOXIDE 

Two  basic  criteria  have  been  applied  in  estimating  the  fugacity  of 
carbon  dioxide  in  the  mineralizing  systems,  namely:  (a)  if  the  fugacity 
of  carbon  dioxide  is  more  than  80  to  100  atmospheres,  almost  all  the 
fluid  inclusions  would  have  shown  C0^  phase  volumes  as  well  as  gas 
inclusions.  Such  features  are  notably  absent  in  those  sections  examined, 
(b)  in  such  a  case,  if  the  fugacity  of  CO^  is  below  50  atmospheres  the 
effect  would  be  negligible  in  computations  of  the  molalities  of  the 
ionic  species,  as  log  values  would  be  mostly  employed.  However, 
calculations  were  made  to  observe  the  variations  in  the  computed  results, 
but  the  effect  is  almost  negligible. 


5.9  MOLALITIES  OF  SULFUR  SPECIES  IN  AQUEOUS  SOLUTIONS 

After  calculating  the  fugacities  of  sulfur,  oxygen,  water  vapor  and 

CO^,  the  fugacities  of  other  gases  such  as  ,  CO,  CH^  and  H^S  were 

calculated  by  applying  Holland's  data,  relative  to  the  major  species. 

According  to  Helgeson  (1969),  the  major  aqueous  species  in  acid,  chloride- 

rich  solutions,  include  HC0  ,  HS0  ,  HS,  H  CO  and  CO  H  S  , 

3  4  2  3  2(aq)  2  (aqj 

KC1,  NaCl,  HC1,  CaSO, ,  NaSO, ,  H.SiO.  and  chloride  and  sulfate  complexes 

4  4  4  4 

of  the  ore-forming  metals.  According  to  Quist  et  al.  (1963)  bicarbonate 
and  bisulfate  complexes  of  sodium  and  potassium  are  very  insignificant 
below  300°C.  Owing  to  the  fact  that  the  necessary  thermodynamic  data, 
as  well  as  experimental  data,  is  not  presently  available,  certain 
assumptions  have  to  be  taken  into  consideration  for  computational 
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convenience,  the  validity  of  which  will  be  discussed  in  a  subsequent 
section.  The  computed  results  of  the  molalities  of  the  major  species 
for  varying  physical  parameters  are  presented  in  Tables  5  and  6. 


5.10  pH  OF  THE  SOLUTIONS 

High  temperatures-pH  controls-CO^  equilibria  and  mica-quartz- 
kaolinite  buffer  relationships  would  enable  the  construction  of  diagrams 
for  a  given  ionic  strength  of  hydrothermal  temperature  at  any  specified 
temperature.  Firstly  the  carbonate  buffer  could  be  operated  on  the 
following  equations: 


nn  (app)  =  B.f 


55.5 


•2C03 


C02  \ 


C0„  C°2 

z 


(2)  H^  =  Henry's 
law 


constant 


H„C0 


2  3  ^ 


HCO. 


H  +  HC£L  K. 

J  i  1 


H+  4-  C03  K2 


aTT+  aTT„ 

H  HC  0„ 


Vs 


aH+  aC0. 


^CO. 


(3) 


(4) 


(CaC03) 


Kc  =  aCa++  aC0. 


(5) 


h2o  — 


H+  +  OH' 


aH+  ’  a0H 


(6) 


211. 


Ca 


0. 


™(CaV 


K^Kc  .  B  .  f  (C02> 


4K  .  (  f  )  .  r 
4  2  V  HC02;  Ca 
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Substituting  and  dividing;  from  the  above  relationships: 


(aH+)3 


(9) 


Relationship  between  pH  and  f(CC>2)  could  thus  be  calculated  from 
Eqn  8.  The  activity  coefficients  could  be  calculated  from  the  relation¬ 
ship  of  the  so-called  'Delta  approximation'  (Helgeson  1969).  The  Delta 
approximation  is  based  upon  the  assumption  that  the  effect  of  ion 
hydration  upon  individual  ion  activity  coefficients  calculated  from 
electrostatic  theory  can  be  represented  by  the  activity  coefficient  of 
a  neutral  species  in  the  supporting  electrolyte.  Helgeson's  selection 
of  t^e  CO^  molecule  as  the  hydrated  molecule  H^CO^,  has  been  applied  as 
the  neutral  species.  Thus,  Delta  approximation  is  expressed  in  the 
following  Eq: 


log  ± 


-Zi2  A  (T)  I1/2 


(I,  T) 


1  +  a°  B  (T)  I 


1/2  +  log  r  H^CO 


(10) 


23(1,  T) 


where  Y\  is  the  activity  coefficient  of  i  th  species,  ^  is  the 

activity  coefficient  of  H  CO  in  the  supporting  electrolyte,  a°  is  the 
effective  diameter  of  the  ion  (assumed  to  be  independent  of  temperature) , 
I  is  the  ionic  strength  of  the  solution,  A  and  B  are  the  Debye-Hdckel 


parameters,  expressed  by: 


A  (T) 


B  (T) 


1.8246  x  106  f>1/2  H20  (T) 

3/2 

[e  (T)  T]J/Z 
8  1/2 

50.29  x  10  j)  H20  (T) 

1/2 

[G  (T)]1/z 


(ID 


(12) 


" 
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where  p  is  the  density  cf  the  fluid,  and  0  is  the  dielectric  constant 
of  water.  Values  of  A  and  B  are  taken  from  Helgeson  (1967c).  In  the 
plot  of  Fig  18  the  relationships  between  temperature-pH  and  fCCO^)  have 
been  given.  Secondly  the  silicate  buffer  could  be  expressed  by  two 
equations : 

3  KA1  Si30g  +  2H%^  KAl3Si3010(OH)2  +  2K+  +  6(Si02)qk{_z  (13) 

(k  mica) 

2  KAl3Si301()(0H)2  +  2H+  +  3^0  ^±3Al2Si203 (OH)  +  2K+  (14) 

(k  mica)  (kaol) 


When  K  feldspar  -  k  mica,  k  mica  -  kaolinite  buffer  systems  are 
taken  into  consideration,  the  following  simplified  equilibrium  relations 
could  be  applied: 


V 


[(aK+)2/K]1/2 


TI  0  ^  1.0 


(15) 


By  varying  pH  in  the  above  equation,  a  relationship  could  be  represented 
between  temperature  and  pH.  When  these  plots  are  superimposed  on  the 
carbonate  buffer  plot,  the  composite  plot  is  shown  in  Fig.  18,  and  the 
position  of  the  equilibrium  for  sericite-quartz-carbonate  is  selected 
at  the  given  temperature  which  would  satisfy  the  maximum  number  of 
parameters  in  question.  However,  the  variations  in  ionic  strengths  would 
show  very  minor  shifts  in  the  equilibrium  positions.  From  the  plot 
(Fig.  18)  the  pH  of  the  ore  fluids  at  350°C,  is  fixed  at  6. 0-6. 5.  This 
infers  therefore  a  very  weak  acidic  to  neutral  pH  for  the  ore-bearing 


solutions  at  350°C. 


Temperature 


FIGURE  18 


pH 


controls  -  CO 


2 


equilibria  -  and  K-feldspar 


-  K-mica  -  Kaolinite  -  quartz  buffer  plot. 


K-fetdspar  -K-mica,  K-mica  -Kaolinite,  Carbonate  buffer 


lon/c  Strength  3*0 
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5.11  K/Na  RATIOS:  INTERPRETATIONS 

Na/K  (atomic)  ratios  for  various  groups  of  analyses  of  connate 
brines  and  fluid  inclusion  extracts  have  been  compiled  by  Sawkins 
(1968)  and  the  frequency  distribution  of  these  fluids  varies  quite 
significantly  from  one  type  of  solution  to  the  other  (compiled  in 
Fig.  19).  On  the  basis  of  these  plots  the  nature  of  the  Yellowknife 
ore-bearing  solutions  were  interpreted.  The  distribution  of  Na/K  vs. 
frequency  (Fig.  19)  shows  that  the  Na/K  ratios  in  connate  waters  are 
greater  than  those  of  sea  water,  whereas  the  reverse  is  true  for  fluid 
from  inclusions  and  volcanic  spring  waters.  An  overall  comparison 
of  the  results  of  the  present  study  with  the  published  frequency 
diagrams  reveals  a  similarity  between  the  providencia  fluid  inclusion 
determinations  on  quartz  and  the  present  determinations  on  quartz  from 
the  Campbell  shear  zone. 

The  variations  of  K/Na  ratios  in  natural  fluids,  in  equilibrium 
with  systems  containing  igneous  melts  have  been  given  by  White  (1968) 

(Fig.  20).  Some  of  the  experimental  results  of  Gammon  (1968)  have 
been  superimposed  on  the  plot  of  White  (1968)  at  higher  temperatures. 

This  plot  signifies  the  necessary  conditions  for  the  fluids  at  any 
given  temperature  to  be  in  equilibrium  with  igneous  melts.  Any  point 
which  lies  far  away  from  the  equilibrium  line  (shown  in  White  1968) 
could  be  explained  in  terms  of  environments  and  history  of  ore  deposition. 
In  other  words,  selective  concentration  or  depletion  of  either  K  or 
Na  could  be  explained  on  the  basis  of  the  geological  history  and 


evidences  in  the  region. 
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Two  possible  explanations  could  be  considered  in  the  light  of 
the  existing  experimental  data.  Firstly,  higher  K/Na  ratios  would 
probably  indicate  initially  higher  temperatures  (Rye  and  Haffty 
1969) .  Had  the  hydrothermal  fluids  been  in  contact  with  the  granitic 
rocks,  which  were  at  temperatures  higher  than  the  temperatures  of 
ore  deposition,  existence  of  such  critical  higher  temperatures  could 
be  easily  explained.  Except  for  the  higher  K/Na  ratios  and  a  possible 
source  from  the  western  granodiorite  intrusion  (as  a  geological 
factor)  there  is  no  positive  evidence  to  support  higher  temperatures 
in  the  initial  stages,  i.e.  secondary  inclusions  should  indicate 
fairly  high  temperatures.  Secondly,  if  K  is  introduced  during  a 
process  of  latter  metamorphic  and  metasomatic  effect,  it  is  quite 
likely  that  the  analyses  of  the  secondary  inclusions  would  have 
shown  higher  K  concentrations.  However,  it  is  highly  impossible  at 
this  stage  to  separate  primary  and  secondary  fluid  inclusions  when 
a  sample  of  25  grams  (approximately)  is  prepared.  A  third  factor, 
though  indirect,  would  have  contributed  a  higher  contamination, 
namely  presence  of  fine  particles  of  muscovite  in  the  form  of 
dust,  associated  with  quartz.  This  could  not  have  been  avoided 
under  any  circumstances.  Higher  K/Na  ratios  in  the  fluids  could 
be  expected  when  a  salinity  of  the  initial  fluids  in  higher  and 


FIGURE  19 


Compiled  frequency  diagrams  of  Na/K  ratio  distributions  in 
different  natural  aqueous  solutions. 

a)  Na/K  ratios  of  connate  waters  and  field  brines 

b)  Na/K  ratios  of  subsurface  waters  of  the  Illinois  Basin 

c)  Na/K  ratios  of  subsurface  water  from  the  Michigan  Basin 

d)  Na/K  ratuios  from  inclusion  fluids  in  Mississippi 
Valley-type  deposits 

e)  Na/K  ratios  of  Providencia  inclusions 

f)  Na/K  ratios  of  spring  waters  from  volcanic  environments 

g)  Na/K  ratios  of  Con  Mine  inclusion  fluids  from  quartz 
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Fig  19 a  (Top).  Na/K  ratios  of  connate  waters  a*'1'; 
field  Sics.  Waited  »vcn,ae  =  «.3.  Data  Iron. i  " 
Hem  ami  Waiinij  lalili-s  12,  13  and  1-1  (1163).  (  "■' 

with  less  than  5,000  ppm  total  ionic  concentration 


Fie.  13c.  (Bottom), 
from  the  Michigan  B 
from  Graf  ct  al.  (1966). 


Xa/K  ratios  of  suhslirtace  ’ 
Basin.  Weighted  average  ~  94.6. 
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Fir..l9d  (Left).  Na/K  ratios  from  inclusion  fluids  in 
Msvissipi’i  Valley- type  deposits.  Average  (not  weighted*) 
?(>.9.  Data  fiom  I  Jail  and  1'iiedman  (1964)  and  Saw  Kins 
l‘i(/ia). 

F  10.13 c  (Middle).  Xa/K  ratios  of  I’rovidencia  inclu-ions. 
.lighted  average  Na/K  :  4.4.  Data  from  Rye  and  llalit) 
:96s). 

Fin.  19 f  (Right),  Na/K  ratios  of  spring  waters  from 
Volcanic  envii onirn  nts.  Weighted  average  14.1.  Data  from 


White,  Hem  and  Waring,  tables  17,  IS,  19,  70  and  21  (1963). 

Note:  Many  of  these  waters  exhibit  low  total  iuiiic  eon 
eenlialion  in  comparison  with  sidiimnlaiy  brint  s  and  Itnid 
inclusions.  This  dilution  is  the  icsult  of  admixtme  witli 
groundwater  of  very  low  ionic  concentration.  I  In  Na/K 
ratios  are  thus  considered  to  he  reasonably  representative 
of  waters  containing  higher  ionic  concentrations  in  the 
deeper  parts  of  volcanic  spting  systems. 


FIGURE  2  0 


Ratios  of  K/Na  in  fluid  phase  vs  temperature,  from 
experimental  data  and  natural  waters  [modified  from 


White  (1968)  and  Gammon  et  al.  (1968)] 


GAMMON  et  al  (1968) 

Chary  (Granite  melts,  1.4  -  2.4  Kb) 

(present  study) 
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subsequent  reactions  might  change  the  salinities  locally  during  the 
process  of  deposition.  Higher  K/Na  ratios  indicate  that  the  hydro- 
thermal  solutions  have  not  mingled  with  sea  water  or  any  other  percolating 
solutions  during  their  transport  through  the  Shear  Zone.  Otherwise  the 
K/Na  ratios  would  have  been  much  lower  than  they  are  at  present.  Thus 
a  more  probable  interpretation  of  the  K/Na  results  would  be  that  the 
chemistry  of  the  solutions  did  not  change  during  the  transport  of  the 
ore-bearing  fluids  through  the  Shear  Zone;  i.e.,  ore  deposition  has  taken 
place  by  magmatic  or  juvenile  waters  and  little  contamination  has  occurred. 
Consistency  in  the  K/Na  ratios  in  all  the  samples  might  indicate  that 
the  temperature  did  not  decrease  very  much  during  the  depositional  history. 

5.12  ISOTOPIC  COMPOSITION  OF  SULFUR  SPECIES 

The  isotopic  composition  of  sulfides  deposited  from  aqueous  solutions 
is  greatly  influenced  by  the  dominant  sulfur  species  present  in  the 
aqueous  solution  at  the  given  condition.  The  chief  sulfur  species  in 
aqueous  solutions  at  temperatures  of  300-360°C  are  H^S,  HS  ,  S  ,  So^, 

HSO. ,  KSO, ,  NaSO. ,  KHSO.,  NaHSO.  and  HnS0.  are  minor  (Ohmoto,  in  press). 

The  isotopic  composition  of  total  sulfur  in  the  fluids  with  respect  to 
the  S  from  the  theoretical  considerations  of  Sakai  (1968)  may  be  expressed 
as  follows: 


S  s3A 


£  S 


(SS3^H2S  .  XH2s)  +  (£S3\s  .  XHS")  +  (SS34=  .  xg=)  + 

s 


(SSSOy  *  V  +  (SSHSOt  *  ^SO--*  +  (§SKSO-  *  XKSO-) 
4  4  4  4  4  4 


34 

KHSO 


34 

NaHSO 


(16) 
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where 

and  X. 
1 

sulfur 


s  34 

d S  is  the  isotopic  composition  of  i  species  in  aqueous  solution, 
is  the  mole  fraction  of  the  species  i  with  respect  to  the  total 
species . 


X. 
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m 


i/("H2S 


"HS 


-  +  V  +  “so; 
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+ 


“hso;  + 

4 


(17) 


The  isotopic  compositions  of  H0S,  HS  ,  and  SO.  were  theoretically 

z  4 

calculated  (Sakai  1968)  and  related  to  S  in  terms  of  the  reducing 

function  ratios.  At  present,  since  the  activity  coefficients  and 

isotope  compositions  of  these  minor  species  are  not  available,  the 

published  and  calculated  values  of  HSO^  and  SO^  have  been  employed  in 

the  computations.  Thode  and  Monster  (1965)  suggested  that  there  is  only 

a  very  small  fractionation  between  SO.  and  sulfate  minerals.  Ohmoto 

4 

(personal  Communication)  has  assumed  the  isotopic  compositions  of  HSO^ , 

KSO  ,  NaSO. ,  KHSO. ,  NaHSO.  to  be  of  the  same  order  as  SO. . 

4  4  4  4  4 

Sufficient  research  has  been  done  to  demonstrate  that  chemical 
equilibrium  is  possible  among  coexisting  sulfur  species  and  precipitating 
sulfides  at  any  given  temperature  (Rye  1969;  Kajiwara  1969;  Grootenboer 
1969).  However,  many  laboratory-controlled  systems  remain  to  be 
analyzed  in  order  to  explain  certain  natural  environments  of  sulfide- 
sulfate  precipitation  and  their  mutual  influence  upon  the  isotopic  compo¬ 
sitions  of  the  precipitating  phases.  However,  the  present  existing  iso¬ 
topic  work  on  sulfur  isotope  fractionation  does  enable  one  to  appreciate 
the  natural  phenomena. 

Mole  fractions,  mentioned  above,  may  be  calculated  if  the  salt 
concentration  is  known,  by  employing  the  following  relationships  at  any 
particular  temperature  (Helgeson  1969) . 
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H2s  ^  H+  +  HS 

Ki 

=  V  •  aHS" 

aH2s 

HS"  ^  H+  +  S= 

K2 

=  V  •  V 

aHS" 

S“  +  2(02)  —  SO^ 

K3 

(asoT) 

=  4 

V  '  f(02)2 

hso7  h+  +  so7 

4  4 

K4 

a  +  •  asoT 

=  4 

aHSo7 

4 

kso7  ^  k+  +  so7 

4  4 

K5 

aK+  *  aso7 

=  4 

aKS07 

4 

NaSOT  ^  Na+  +  SOT 

4  4 

K6 

a  +  .  a  = 
Na  SO, 

=  4 

^aSOj 

4 

KHSO,-^  K  +  HSO, 

4  4 

K7 

aK+  *  aHS07 
=  4 

aKHS0. 

(18) 
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(20) 
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NaHSO.  ==iNa+  +  HSO. 
4  4 
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aNa  *  aHSO, 


^aHSO, 
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H+  +  OH 
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H  SO. 
2  4 


K 


4 


4 


10 


(27) 
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a . 
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where  K^,  K^,  ...  are  the  equilibrium  constants  for  reactions  at  the 
given  temperature  (Helgeson  1969)  and  a^,  ^ ,  nu  represent  the  activity, 


activity  coefficient  and  molality  respectively.  Computational  data 
in  the  calculation  of  the  activity  coefficients  of  different  species,  by 
using  Delta  approximation  have  been  given  in  the  appendix  2. 

Recent  work  on  the  isotopic  composition  of  sulfur  species  (Ohmoto, 
personal  communication)  indicates  that  the  factors  influencing  the 
isotopic  composition  of  various  species  and  mineral  phases  being  precipi¬ 
tated  from  these  solutions,  are  temperature,  fugacity  of  oxygen,  hydrogen 
ion  concentration,  ionic  strength  of  the  solutions,  concentration  of  K 
and  Na  in  solution,  and  isotopic  composition  of  total  sulfur. 

For  the  present  study  on  the  hydrothermal  model,  with  the  established 
data  derived  from  the  experimental  work,  theoretical  calculations  and 
other  reasonable  assumptions,  the  effect  of  these  parameters  on  the 
molalities  of  the  species  and  on  the  isotopic  composition  of  the  sulfides 
in  aqueous  solutions  have  been  examined.  The  molalities  of  the  sulfur 


given  in  the  tables  6  and  7  for  the  parameters  specified.  The  major 
computational  procedures  to  calculate  the  individual  molalities,  thereby 
calculating  the  isotopic  composition  of  sulfides  in  aqueous  solution  are 
mainly  based  on  equations  16  to  28.  The  computational  procedures  are 


outlined  in  Appendix  3. 

From  Tables  6  and  7,  it  can  be  shown  that  a  decrease  in  temperature 
would  move  the  stability  fields  of  minerals  towards  lower  fugacities  of 
oxygen  and  sulfur.  The  molalities  of  the  constituents  would  also  vary 
considerably. 

Within  the  range  of  interest,  the  molalities  of  H^S,  HS  increase 
slightly,  whereas  the  molalities  of  the  rest  of  the  sulfur  species 
decrease  considerably  with  a  lowering  of  the  temperature  and  the  retention 
of  remaining  factors  at  a  constant  level,  except  a  corresponding  change 
in  the  fugacities  of  oxygen  and  sulfur  with  reference  to  the  stability 
fields  of  sulfides. 

The  molalities  of  the  sulfur  species  are  not  directly  related  to 
the  pressure  except  that  the  fugacity  of  water  vapor  in  the  system  is 
related  to  the  total  pressure  by  a  fugacity  coefficient.  However,  the 
effect  of  the  partial  pressure  of  water  vapor  is  quite  small  in  computing 
the  molalities  of  the  constituents.  Thus  all  the  computations  are  for 
an  average  pressure  of  500  atmospheres. 

At  constant  temperature,  ionic  strength,  hydrogen  ion  concentration, 
and  pressure,  decreasing  both  fugacities  of  oxygen  and  sulfur  would 
result  in  a  decrease  of  the  molalities  of  H^S,  HS  and  S'and  other  sulfate 
species.  But  at  constant  factor  of  f(C>2),  a  decrease  in  f(S2>  would 
decrease  the  molalities  of  H2S,  HS*,  S'and  other  sulfate  species  and  at 
constant  f(S2),  a  decrease  in  f(02)  would  increase  the  molalities  of 

other  sulfide  species. 

Ionic  strength  has  no  considerable  effect  upon  the  molalities  of 
the  constituents.  Changes  in  the  molalities  are  so  small  that  the  indi¬ 
vidual  percentages  of  the  species  would  not  vary  considerably.  Since 
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TABLE  6  and  TABLE  7 


Computed  molalities  of  the  sulfur  species 
in  aqueous  solutions  at  350°C  and  300°C. 
6.0  (rest  of  the  parameters  specified). 


considered 
pH  3.0  and 


N.B.:  A  pH  of  3  is  selected  to  show  the  variations  in 
the  molalities  of  the  species  with  pH 


Molalities 
log  values 

I  =  3.0 
pH  =  3.0 

log  f(S2) 

in 

350 

°C 

300 

°C 

log  f(02 

-10.5 

)=-31 

-8.5 

log  f(o2: 

-10.5 

>=-29 

-8.5 

log  f (0 

-11.5 

2)-34 

-9.5 

log  f(02 

-11.5 

)=-32 

-9.5 

■k 

h2s 

-3.22 

-2.22 

-4.22 

-3.22 

-2.76 

-1.76 

-3.76 

-2.76 

HS 

-7.71 

-6.71 

-8.7. 

-7.71 

-7.14 

-6.14 

-8.14 

-7.14 

S 

-9.24 

-8.24 

-10.24 

-9.24 

-10.51 

-9.51 

-11.51 

-10.51 

SO, 

4 

-12.19 

-11.19 

-9.19 

-8.19 

-13.88 

-12.88 

-10.88 

-9.88 

HSO, 

4 

-9.49 

-8.49 

-6.49 

-5.49 

-11.60 

-10.60 

-8.60 

-7.60 

KSO, 

4 

-11.37 

-10.37 

-8.37 

-7.37 

-13.38 

-12.38 

-10.38 

-9.38 

NaSO, 

4 

-10.99 

-9.99 

-7.99 

-6.99 

-13.07 

-12.07 

-10.07 

-9.07 

KHSO, 

4 

-8.67 

-7.67 

-5.67 

-4.67 

-14.67 

-10.61 

-8.61 

-7.61 

NaHSO, 

4 

-8.29 

-7.29 

-5.29 

-4.29 

-11.29 

-10.29 

-8.29 

-7.29 

H  SO, 

-13.93 

-12.93 

-10.93 

-9.93 

-15  72 

-14.72 

-12.72 

-11.72 

2  4 


*  log  values 


NOTE:  All  these  values  are  relative  quantities,  not  absolute  values 


I  =  3.0 

350° 

C 

300 

°C 

pH  =  6.0 

log  f(02) 

=-31 

log  f(o2: 

>=-29 

log  f(02 

)=-31 

log  f(02 

)=-29 

log  f(S2) 

-10.5 

-8.5 

-10.5 

-8.5 

-11.5 

-9.5 

-11.5 

-8.5 

* 

H2s 

-3.22 

-2.22 

-4.22 

-3.22 

-2.76 

-1.76 

-3.76 

-2.76 

HS” 
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-3.71 
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-4.71 

-4.14 

-3.14 

-5.14 

-4.14 

S= 

-3.24 

-2.24 

-4.24 

-3.24 
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-3.51 

-5.51 

-4.51 

soT 

4 

-6.19 

-5.19 

-3.19 

-2.19 

-7.88 

-6.88 

-4.88 

-3.88 

HSO. 

4 

-6.49 

-5.49 

-3.49 

-2.49 

-8.60 

-7.60 

-5.60 

-4.60 

ksoT 

4 
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2  4 

*  Molalities  in  the  log  values 


NOTE:  All  these  values  are  relative  quantities,  not  absolute  values 
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this  effect  would  change  the  percentages  imperceptibly,  for  the  present 
computational  purposes,  an  ionic  strength  of  3.0  is  selected  (Helgeson 
1969) . 

At  constant  temperature,  pressure,  ionic  strength,  fugacities  of 
oxygen  and  sulfur,  a  decrease  in  the  hydrogen  ion  concentration  or  an 
increase  in  pH  would  not  bring  about  a  change  in  the  molalities  of 
within  the  range  of  interest,  i.e.,  3.0  to  6.0,  but  a  considerable 
increase  in  the  molalities  of  HS*,  S  and  all  other  sulfate  ions. 

With  regard  to  the  K/Na  ratios,  different  values  for  these  ratios 
would  influence  the  molalities  of  potassium  and  sodium  complex  sulfates 
and  would  thus  change  their  individual  molalities.  However,  this  effect 
is  quite  negligible  when  considering  all  other  major  factors. 

5.13  THE  INFLUENCE  OF  PHYSICO-CHEMICAL  FACTORS  ON  THE  ISOTOPIC 
COMPOSITION  OF  SULFIDES 

As  discussed  previously,  at  any  given  temperature,  the  hydrogen  ion 
concentration  and  fugacity  of  oxygen  are  the  major  factors  influencing 
the  isotopic  composition  of  sulfides  in  aqueous  solutions  as  the  molalities 
of  the  individual  species  are  greatly  modified  by  changing  the  pH  and 
fugacity  of  oxygen.  According  to  Sakai's  theoretical  calculations  when 
S"  is  assumed  to  be  zero  at  all  instances,  the  isotopic  composition  of 
H^S,  HS  ,  pyrite  and  SO^  would  be  4.8,  3.7,  4.5  and  23.5  permil  heavier 
than  the  isotopic  composition  of  S  at  350°.  Based  upon  this  assumption, 
the  isotopic  compositions  of  sulfides  in  aqueous  solutions  have  been 
calculated  at  varying  pH  and  f(C>2)  and  the  plot  presented  in  Figure  20. 

From  these  calculated  isotopic  compositions  the  isotopic  composition  of 
pyrite  would  have  varied  considerably.  The  ranges  of  oxygen  fugacity  were 
fixed  with  reference  to  Holland's  data.  Actually  the  lower  limit  for  pH 
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is  fixed  from  Helgeson  (1967)  of  the  aqueous  solutions.  The  upper  limit 
is  taken  from  the  silica,  mica  and  carbonate  buffer  (Fig.  18).  The  plot 
clearly  shows  that  the  isotopic  composition  of  pyrite  would  assume  posi¬ 
tive  values  towards  the  end  phases  of  deposition  if  the  physico-chemical 
conditions  varied  accordingly.  It  is  not  very  essential  that  the  temp¬ 
erature  vary  considerably  to  bring  about  major  changes  in  the  isotopic 
composition  of  the  species.  Even  slight  changes  would  bring  considerable 
changes  in  the  isotopic  composition. 

There  are  thus  three  factors  that  clearly  indicate  that  the  initial 
pH  would  have  been  greater  than  3.0.  As  Helgeson  stated  (1967)  that  a  pH 
of  3.0  is  unrealistic  geologically  if  carbonate  minerals  are  considered 
to  have  been  deposited  simultaneously  with  the  gold  quartz  bodies. 

Secondly,  if  the  pH  is  considered  to  be  3.0  at  the  specified  maximum  of 
fugacity  of  oxygen,  the  fugacity  of  sulfur  is  comparatively  very  high. 

If  the  corresponding  fugacities  of  oxygen  and  sulfur  are  taken  into 
consideration,  the  pH  of  the  solutions  would  have  been  around  6.0  and  the 
isotopic  composition  of  aqueous  sulfides  would  be  -23.5%»  Had  pyrite  been 
precipitated  from  these  aqueous  solutions  of  isotopic  composition  -23.5%, 
it  would  show  an  isotopic  composition  of  -19.0  at  a  fugacity  of  oxygen  of 
10  exp.  -29.0  and  f(S  )  10  exp.  -8.5.  Even  a  slight  decrease  in  fugacities 
of  sulfur,  oxygen  or  pH  would  shift  the  isotopic  composition  of  aqueous 
sulfides  to  much  heavier  values  and  consequently  resulting  in  much  heavier 
isotopic  composition  for  pyrite. 

In  this  critical  region,  increase  of  pH  by  1  unit,  decrease  of  the 
fugacity  of  oxygen  by  two  units  (log)  and  the  fugacity  of  sulfur  by  two 
log  units  would  change  the  isotopic  composition  of  sulfides  in  aqueous 
solution  by  23.0%.  and  a  corresponding  shift  could  be  expected  in  the 


FIGURE  21 


Computed  plot  showing  the  relationship  between  the 

34 

fugacity  of  oxygen  and  S  %0 on  an  iso-pH  projection 
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isotopic  composition  of  pyrite  and  other  sulfides  deposited  from  the 
solutions  (Fig.  21).  If  galena  and  sphalerite  are  assumed  to  have  been 
co-precipitated  along  with  pyrite,  the  corresponding  minimum  and  maximum 
isotopic  composition  for  galena  would  be  as  follows: 

&  s34  %. 


s 

Pyrite 

Sphalerite 

Galena 

£H 

-21.0 

-16.5 

-18.0 

-20.2 

3.0 

-0.3 

+4.2 

+2.7 

+0.5 

6.0 

This  would  have  been  a  more  probable  reason  for  variations  in  the 
isotopic  compositions  of  different  sulfides  collected  from  the  Campbell 
Shear  Zone  system.  Temperature  would  not  have  been  a  major  factor  as 
the  variations  in  the  isotopic  compositions  are  relatively  small.  From 
the  experimental  data,  it  could  be  interpreted  that  fugacity  of  oxygen 
would  not  have  changed  quite  noticeably,  otherwise  an  abrupt  change  in 
the  isotopic  composition  would  have  been  occurred  from  a  very  negative 
permil  value  to  more  positive  values  (Fig.  19) . 

5.14  OXYGEN  ISOTOPES:  INTERPRETATIONS 

From  the  oxygen  isotopic  compositions  determined  on  CO^  gas 
extracted  from  the  carbonates,  the  oxygen  isotopic  compositions  of  the 
quartz  and  the  mineralizing  waters  were  calculated,  assuming  equilibrium 
among  the  carbonates,  quartz  and  ore-fluids  at  a  given  temperature. 
Samples  with  high  and  low  values  have  been  selected  for  calculating 
extreme  values.  Clayton's  (1967)  empirical  equations  have  been  employed: 
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‘(quartz-water) 


=  1000  In  =  3.57  (106  T  2)  -2.73 


6  -2. 


A(calcite-water)  =  1000  ln  *  =  2*80  <10  T  >  "3 


1000  In °c  .  <■  18  Cn18 

min-water  =50.  -  50 

mm  water 


.40 


(29) 

(30) 

(31) 


The  values  of  hypothetical  H^O  coexisting  with  a  given  mineral 
assemblage  (Viz.  quartz  and  calcite)  at  a  given  temperature  were  calcu¬ 
lated,  based  upon  the  mineral-water  fractionation  curves. 
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'quartz 


8.7-10.' 


5. 0-6. 8%.  21.4-19.6%, 


temp 


350°C 


As  no  temperature  gradient  during  ore  deposition  has  been  recognised, 

it  might  be  assumed  that  there  were  no  changes  in  the  isotopic  composition 

during  cooling.  In  the  observed  mineral  assemblages,  quartz  is  invariably 
18 

rich  in  0  because  of  isotopic-exchange  reactions,  where  the  frequency 

18  16 

shifts  on  isotopic  substitution  of  0  for  0  are  similar  to  two  compounds, 
18 

0  should  be  concentrated  by  the  compound  whose  oxygen  atoms  have  the 


J8 


highest  vibrational  frequencies. 

Clayton  and  Epstein  (1961)  found  that  calculated  £  O"1^  values  for 
hydrothermal  waters  have  a  mean  of  6%,  O'Neil  and  Clayton  (1964)  indi¬ 
cated  that  the  isotopic  composition  of  vapor  phases,  presumed  to  have 
been  in  equilibrium  with  various  igneous  rock  types,  falls  in  a  narrow 
range  near  8%»  The  Providencia  fluids  (Rye  1966)  have  been  reported  to 
be  magmatic  waters.  The  isotopic  composition  obtained  for  the  Gilman 
ore  body  waters  was  8 . 6-6 . 8%,  which  were  interpreted  to  be  normal  magmatic 
waters  (Engel  et  al.  1958).  An  analogy  could  be  drawn  between  the  Gilman 
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ore  body  and  the  Campbell  shear  system  considering  the  geological  setting 
of  the  deposits.  Furthermore  the  following  generalizations  were  inter¬ 
preted  based  on  the  experimental  data  from  a  number  of  ore  deposits. 
Taylor  (1967)  (1)  magmatic  waters  are  relatively  uniform  in  0^/0^ 
ratio,  with  values  of  7.5  to  9.0  permil.  (2)  Hydrothermal  waters 
very  commonly  have  S  values  similar  to  those  of  magmatic  waters  with  a 

very  few  exemptions.  (3)  Metamorphic  waters,  particularly  those  in 

18 

contact  with  meta  sediments,  generally  are  higher  in  0  than  magmatic 
waters  (10-14%).  (4)  Oceanic  water  is  very  uniform  at  -0.5  to  +0.5^0 

(5)  Meteoric  and  surface  fresh  waters  generally  show  -10  to  -4%„ 


5.15  CARBON  ISOTOPES:  INTERPRETATIONS 


Assuming  that  the  major  fluid  components  carrying  carbon  are  CO^ 
and  CH^  in  aqueous  solutions,  the  fugacities  of  CO^  and  CH^  have  been 
utilized  in  computing  the  total  carbon  isotopic  composition  by  using  the 
following  relationship: 
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The  parameters  used  in  the  calculations  are  the  fugacity  of  carbondioxide , 
the  fugacity  of  methane  (in  atmospheres),  the  total  pressure,  the  temp¬ 
erature  and  the  derived  isotopic  composition  of  carbondioxide  from 
carbonates.  The  fugacity  of  CH^  may  be  calculated  from  the  relationship: 


2H^0  +  (C02)  ^  CH4  +  2(02)  [Holland  1965] 

The  isotopic  compositions  of  C02  and  CH^ ,  relative  to  CaCO^  are  taken 

13 

from  Bottinga's  (1969)  calculated  values  for  1000  In  c<  for  C  exchange 
relationships  among  carbondioxide,  calcrte  graphite  and  methane. 

The  computed  values  of  isotopic  compositions  of  total  carbon  in 
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aqueous  solution,  range  from  -4.4  to  -0.6^o  for  two  extreme  values 
obtained  on  the  carbonates,  viz.  (-4.2 %0  and  -0.5%o).  With  the  limited 
amount  of  data,  it  is  not  possible  to  establish  any  trends  of  &  C  during 
calcite  precipitation.  The  most  obvious  temperature  dependent  reaction 
during  calcite  deposition  was  the  isotopic  fractionation  of  carbon 
between  the  hydrothermal  fluids  and  precipitating  calcites.  If  the 
calculated  values  are  correct ,  they  fall  near  the  range  recorded  for 

materials  believed  to  represent  juvenile  carbon.  Craig  (1953)  estimated 

13 

a  b  C  of  terrestrial  carbon  to  be  -7%  from  material-balance  consider- 

13 

ations  of  the  carbon  cycle  in  nature.  S>  C  values  reported  for  diamonds 

range  from  about  -2.0  to  -9.0%o  and  average  about  -6 %Q  (Wickman  1956). 

13 

8  C  determinations  in  carbonatite  ejecta  in  pleistocene  tuffs  from  the 
Laacher  See  district,  Germany,  which  are  believed  to  represent  primary 
igneous  carbonatite  derived  from  a  gabbroic  magma,  give  a  range  of  -6.6 
to  -8.4%0  (Taylor  et  al.  1965).  Similarly,  Providencia  hydrothermal 
fluids  have  been  interpreted  to  contain  deep-seated  carbon  in  the  fluids 
(Rye  1966). 

The  possibility  of  a  juvenile  source  for  the  carbon  in  the  hydro- 
thermal  fluids  at  Yellowknife  is  in  substantial  agreement  with  the  geo¬ 
logical  picture  of  the  area  having  the  western  granodiorite  intrusive 
as  the  major  source.  The  uncertainties  in  determining  the  fugacities  of 
oxygen  and  water  vapor  or  the  experimental  uncertainties  in  determining 
the  carbon  isotopic  composition  would  be  the  major  factors  of  error. 
However,  detailed  work  is  planned  for  further  studies.  Since  carbon 
isotopic  composition  also  falls  within  the  range  of  magmatic  carbon 
isotopic  composition,  it  is  more  justifiable  to  state  that  the  source  of 
carbon  as  well  as  the  ore  solutions  is  from  a  magmatic  entity . 
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5*16  THEORETICAL  ASSUMPTIONS,  UNCERTAINTIES  AND  SOURCES  OF  ERRORS 

Before  evaluating  the  validity  of  the  results,  it  is  essential  to 
discuss  the  sources  of  errors  in  the  computations. 

1.  In  the  theoretical  calculations  of  Sakai  (1968)  the  equilibrium 
among  the  given  set  of  sulfur  species  is  assumed,  but  in  nature,  reactions 
tend  to  reach  equilibrium. 

2.  The  computed  values  of  the  isotopic  compositions  of  the  sulfur 
compounds  (log  f)  are  approximations. 

3.  In  nature,  closed  system  reactions  are  rare. 

4.  Since  the  isotopic  compositions  of  HSO^ ,  KSO^,  NaSO^,  KHSO^, 
NaHSO^  at  a  given  temperature  are  not  theoretically  calculated,  for  the 
present  purposes  it  is  assumed  that  they  are  almost  of  the  same  order 
of  SO^  (Ohmoto,  personal  communication). 

5.  Since  bicarbonate  and  bisulfate  complexes  of  sodium  and  potassium 
are  very  insignificant,  they  are  not  considered  in  the  computations. 

6.  The  molalities  of  the  different  species  calculated  are  relative 
quantities  rather  than  absolute  values. 

7.  Equilibrium  reactions  are  assumed  among  hydrothermal  fluids, 
quartz,  and  calcite. 

8.  The  major  carbon  species  considered  are  CO^  and  CH^  in  equi¬ 
librium  within  the  system. 

9.  Uncertainties  in  the  equilibrium  constants  and  activity  co¬ 
efficients  . 

5.17  SOURCE  OF  ORE  CONSTITUENTS 

According  to  White  (1968)  the  following  sources  for  a  given  ore 
deposit  may  be  discussed  in  the  light  of  the  geology  and  geochemistry  of 
the  area:  (a)  metals  dispersed  in  products  of  weathering  and  sediments, 
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(b)  metals  undergoing  metamorphism,  (c)  metals  in  rocks  of  the  deep 
crust  or  the  mantle,  or  (d)  metals  present  in  earlier-formed  ore 
deposits.  Two  of  the  sources  mentioned  above  seem  pertinent  to  the  Con 
Mine  deposits,  namely,  (1)  rocks  being  metamorphosed,  or  (2)  rocks  of 
the  deep  crust  or  mantle.  The  first  case  is  supported  by  Boyle  (1961) 
whereas  Ames  (1962)  ,  Coleman  (1957)  and  McConnell  (1964)  supported  the 
second  theory. 

The  observation  by  McConnell  (1964) ,  of  narrow  veinlets  within  the 
granitic  rocks  West  of  the  Yellowknife  greenstrone  belt,  containing  small 
amounts  of  sulfides,  (pyrite,  minor  galena  and  sphalerite)  and  erratic 
gold  values  is  perhaps  pertinent  to  the  present  discussion.  Green's 
work  (1968)  strongly  supports  the  theory  that  the  Western  Granodiorite 
might  not  have  formed  by  the  granitization  process.  In  the  light  of 
Boyle's  (1961)  lateral  secretion  theory,  it  is  difficult  to  explain  the 
presence  of  the  ore  deposit  in  the  granite  body.  Thus  McConnell  strongly 
feels  that  the  gold-quartz  deposits  in  the  Yellowknife  area  were  formed 
and  localised  due  to  the  ascending  hydrothermal  solutions  through  the 
Shear  Zones,  resulting  in  alteration  haloes  surrounding  the  lenses. 

5.18  THE  CHEMICAL  NATURE  OF  THE  ORE-FLUIDS 

The  chemistry  of  the  ore-fluids  depends  principally  upon  the 
dominent  factor  influencing  the  chemical  composition  as  well  as  the 
dominent  chemical  constituent  of  the  solutions.  Sulfide  complexes  of 
As,  Sb  and  Hg  are  highly  soluble  at  high  pH's  and  at  concentrations  of 
total  sulfide  (Dickson  1961).  The  bisulfate  complexes  of  Cu,  Zn,  Cd, 

Pb,  Ag,  Hg,  Sb  and  As  seem  to  be  adequately  stable  (Barnes  and  Czamanski 
1967)  at  reasonable  pH.  Regarding  the  character  of  the  ore-bearing 
solutions,  Ridland  (1941)  stated  that  the  main  constituents  in  the  ore 
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fluids  are  Ta,  Mn,  Fe,  As,  Sb,  Zn,  Ca,  Cu,  Te,  Au,  Ag,  Fb,  Si,  S  and 
^0*  As  has  been  calculated  from  the  silicate  and  carbonate  buffer 
system  (Fig.  18)  a  pH  of  5. 5-6. 5  would  quite  adequately  explain  the 
isotopic  variations  in  sulfides.  These  weak  acid  solutions  would  have 
influenced  the  formation  of  the  alteration  zones  surrounding  the  lenses. 
Geological  and  isotopic  evidence  shows  that  the  source  was  magmatic  but 
to  obtain  possible  evidence  regarding  this  conclusion  and  contrary  to 
Boyle's  theory  of  lateral  secretion  (1961),  Deuterium,  Sr,  Pb  isotope 
work  is  obviously  necessary  to  facilitate  comparison  between  the  possible 
'parental'  granite  body  and  the  ore  deposits.  Oxygen  isotope  work  on 
carbonates  and  other  silicate  phases  would  hopefully  augment  such  infor¬ 
mation. 

The  solutions  are  thought  to  have  percolated  through  the  Shear  Zone 
and  filled  the  dilation  sectors.  The  mechanism  of  transport  could  be 
expected  as  due  to  (1)  the  regional  geothermal  gradient,  or  (2)  the 
pressure  potential  within  the  ore  solutions. 

Age  of  Mineralization 

The  problem  of  age  of  mineralization  presently  rests  in  the  state  of 
a  logical  speculation.  Ridland  (1941)  related  the  gold  mineralization  in 
the  Yellowknife  area  to  one  of  the  following  possibilities:  (a)  as  the 
gold  deposits  are  found  in  the  youngest  pre— Cambrian  rocks  of  the  district, 
they  may  have  been  formed  during  the  late  pre— Cambrian  orogenies  in  which 
the  younger  sedimentary  rocks  were  folded  and  metamorphosed;  (b)  a  giant 
quartz  vein,  185  miles  Northwest  of  Yellowknife,  containing  pitchblende, 
has  been  dated  by  lead-uranium  ratios  as  late  Ordovician  and  the  mineral¬ 
ization  may  be  related  to  this  period,  and  (c)  the  Franklin  Mountains,  140 
miles  West  of  Yellowknife,  were  formed  during  Laramide  revolution.  It  is 
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possible  that  the  Yellowknife  gold-quartz  deposits  were  formed  during  this 
orogeny.  Robertson  and  Cumming  (1966)  obtained  the  lead  model  date  2900  m.y. 
on  sulfides  from  veins  in  the  metavolcanics ,  and  Green  (1968)  obtained 
K-Ar  and  R-Sr  data  for  the  granodiorites ,  the  date  2610  +  58  m.y.  and 
metavolcanics  2625  ±  60  m.y.  There  is  a  slight  discrepancy  in  the  results 
obtained  by  these  two  approaches.  However,  it  is  suggested  that  the  out¬ 
pouring  of  volcanic  rocks  took  place  only  slightly  before  the  intrusion  of 
granodiorite  batholiths,  to  which  mineralization  in  the  district  could  be 
attributed.  McConnell’s  (1964)  observations  also  support  this  conclusion. 
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CHAPTER  VI 

SUMMARY  -  CONCLUSIONS  AND  INDICATED  FURTHER  RESEARCH 

From  the  data  gained  from  the  sulfur  isotope  fractionation  studies 
and  fluid  inclusion  studies  it  is  concluded  that  mineralization  in  the 
gold  quartz  veins  of  Yellowknife  area  took  place  at  a  temperature  of 
350°-300°C.  The  general  nature  of  the  gold-ore-bodies,  the  persistence 
of  ore  with  increasing  depth  and  the  character  of  the  major  mineral 
constituents,  suggest  that  p-T  conditions  during  the  mineralization  were 
possibly  uniform  and  rank  in  medium  to  high  grade.  In  the  light  of  the 
thermodynamic  data  it  is  shown  that  the  isotopic  composition  of  the 
sulfides  is  a  function  of  physico-chemical  conditions  during  precipitation. 
From  the  microscopic  paragenetic  studies,  it  is  postulated  that  mineral¬ 
ization  most  likely  took  place  in  two  stages,  this  fact  being  augmented 
by  the  isotopic  compositions  of  the  sulfides  (the  earlier  pairs 
exhibiting  lighter  isotopic  composition  compared  with  the  later  pairs) . 

Calculated  oxygen  isotopic  compositions  of  the  hypothetical  ore 
fluids,  as  well  as  the  isotopic  compositions  of  the  carbon  in  the  ore 
solutions,  fall  within  the  limits  of  the  isotopic  compositions  of  the 
respective  constituents  of  the  magmatic  fluids.  Furthermore,  the  K/Na 
ratio  frequency  distribution  diagrams  indicate  a  comparison  with  the 
compositions  of  magmatic  fluids. 

The  probable  concentration  of  hydrogen  ions  in  the  aqueous  solutions 
at  350°C  was  determined  by  adopting  the  silica  -  K-mica  -  carbonate 
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buffer  system.  The  pH  of  the  mineralizing  solutions  was  thus  determined 
to  be  in  the  range  of  5. 5-6.0.  The  mineralizing  fluids  of  pH  5. 5-6.0, 
when  passing  through  the  Shear  Zone,  might  have  reacted  with  the  country 
rock,  as  a  result  of  which  sericite-carbonate  schist  has  formed, 
enveloping  each  quartz  lens. 

The  variation  of  the  molalities  of  different  sulfur  species  in  the 
aqueous  mineralizing  solutions  have  been  demonstrated  by  computed  plots 
and  the  conclusion  is  drawn  that  pH,  fugacity  of  oxygen  and  temperature 
of  formation  would  have  influenced  the  minor  variations  in  the  isotopic 
compositions  of  the  sulfides  in  the  aqueous  solutions  to  a  major  extent 
in  the  case  of  Campbell  shear  system  whereas  the  influence  of  ionic 
strength,  K  and  Na  ratios  and  partial  pressure  of  water  vapor  is 
negligible . 

From  the  techniques  used  in  deciphering  the  nature  of  the  ore 
solutions  and  their  chemistry,  it  is  felt  that  further  research  to  clarify 
some  of  the  unsolved  and  critical  questions,  is  quite  essential. 

Particular  mention  may  be  made  of  more  detailed  work  on  the  stable 
isotope  geochemistry  within  the  alteration  zone;  i.e.,  a  study  of  the 
variations  of  sulfur,  oxygen  and  carbon  isotopic  compositions  in  the  ore 
bodies,  the  surrounding  alteration  zones  and  the  country  rocks,  would 
perhaps  enable  one  to  detect  the  source  of  the  ore  constituents  and  their 
ultimate  mode  of  emplacement.  Similarly  the  following  further  research 
is  suggested:  (a)  electron  microprobe  studies  on  element  partition  during 
the  ore  forming  processes,  viz;  K,  Na,  Sr,  etc., (b)  a  study  of  deuterium 
isotope  variations  to  further  investigate  the  nature  of  the  ore  forming 
processes,  (c)  geochemical  studies  on  the  sericitic  micas  from  the 
alteration  zones,  (d)  the  determination  of  chlorine  concentrations  which 
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might  be  a  geochemical  indicator  of  potential  ore  zones,  and  (e)  based 
on  the  empirical  studies  on  the  solubilities  of  gold,  quartz,  pyrite 
and  calcite,  quantitative  studies  of  the  gold,  enriched  in  each  lens 
from  the  known  data  of  the  dimension  of  the  quartz  lenses. 
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APPENDIX  1 

Thermodynamic  Data  (Holland  1965,  1959) 


A 

B 

C 

D 

Ref 

1. 

3/2[Fe]  +(02)  l/2[Fe304] 

-133900 

41.10 

298-833 

1 

2. 

2 [Fe]  +(S2)  2 [FeS] 

-71820 

25.12 

412-1179 

1 

3. 

4[Fe304]+(02)  6[Fe20  ] 

-114700 

67.8 

300-1200 

2 

4. 

Fe203+(S2)+5/2(O2)  2[FeS04] 

-275000 

147.1 

300-1000 

1 

5. 

2[FeS]+(S2)  2[FeS2] 

-71000 

71.9 

600-1000 

2 

6 . 

[FeS2]+2(02)  [FeS04]+l/2(S2) 

-162600 

57.8 

300-1000 

1 

7. 

2 [Zn]+[ S2 ]  2 [ZnS] 

-127300 

-9.2 

73.4 

298-693 

1 

8. 

2[Zn]+(02)  2 [ZnO] 

-168200 

-13.8 

88.3 

298-693 

1 

9. 

[ZnS]+2 (0  )  [ZnSO . ] 

2  4 

-185710 

-7.37 

105.0 

298-1000 

2 

10. 

[Zn0]+l/2(S2)+3/2(02)  [ZnS04] 

-166200 

80.5 

300-1000 

1 

11. 

2[Pb]+S2  PbS 

-72280 

41.3 

600-1200 

2 

12. 

2[Pb]+(02)  2[Pb02] 

-109900 

-16.1 

100.2 

600-1150 

1 

13. 

[PbS]+2(02)  [PbS04] 

-197000 

84.3 

300-1000 

2 

14. 

[Pb0]+l/2(S2)+3/2(02)  [PbS04] 

-182600 

81.7 

300-1000 

1 

15. 

[s±]+o2  [Si02] 

-210600 

-3.0 

52.22 

298-1700 

1 

16. 

[Si]+(S2)  [SiS2] 

-80000 

39.7 

298-1000 

2 

17. 

[si02]+(s2)  [sis2]+(o2) 

130000 

3.0 

-12.5 

298-1000 

2 

18. 

[c]+(02)  [co23 

-94200 

-0.2 

298-2000 

2 

19. 

[CaS04]+(C02)  [CaC03J+l/2(S2) 

+3/2(02) 

+163500 

-46.3 

300-1000 

1 

20. 

[CuFeS2J+S2  [Cu5FeS4]+4[FeS2] 

-50728 

56.95 

525-775 

2 

21. 

(H20)+1/2(S2)  (H2S)+l/2(02) 

37145 

-0.86 

1.61 

298-1750 

2 

22. 

S2+2(02)  2(S02) 

-173240 

34.62 

298-2000 

2 

1  Holland  (1959) 


2  Holland  (1965) 


Molal  Values  of  the  Dehye-TIuelccI  Fleet  rost  n  I  io  I’a  ra  me  t  cr> 
A  and  II  from  25  to  550"  (.  ami  the  Values  of  tin-  Dieiee.trie  ('.oiimI.iiiI 
(Akcrlof  and  Oshry,  1950)  and  Density  (Keenan  and  Kcvcs,  1950) 
of  Water  Fniployed  in  the  Calculations* 


Temperature 

eC 

A 

D  X  10-' 

<njO 

PrtjO 

25 

0.5095 

0.3284 

78.420 

0.99GS 

30 

0.5144 

0.3292 

76.618 

0.9956 

40 

'  0.5244 

0.3310 

73.125 

0.9919 

50 

0.5354 

0.3329 

69.800 

0.9SS2 

00 

0.5471 

0.3347 

66.629 

0.9833 

70 

0.5596 

0.3366 

63 . 604 

0.9779 

75 

0.5661 

0.3376 

62.144 

0.9750 

80 

0.5729 

0.3386 

60.718 

0.9720 

90 

0.5871 

0.3406 

57.965 

0.9656 

100 

0.6019 

0.3425 

55-.  336 

0.95S0 

110 

0.61S0 

0.3447 

52.826 

0.9512 

120 

0.6347 

0.3468 

50.427 

0.942S 

125 

0.6435 

0.3479 

49.268 

0.9390 

130 

0.6525 

0.3490 

48.134 

0.9346 

140 

0.6715 

0.3512 

45.939 

0.9259 

150 

0.6915 

0.3536 

43.838 

0.9169 

160 

0.7129 

0.3559 

41.824 

0.9076 

170 

0.7354 

0.3583 

39.891 

0.8974 

17^ 

0.7472 

0.3595 

38.953 

0.8921 

180 

0.7595 

0.3608 

38.033 

0.8S70 

190 

0.7851 

0.3632 

36.244 

0. 875S 

200 

0.8127 

0.3659 

34.519 

0.8645 

210 

0.8427 

0.36S6 

32.852 

0.8530 

220 

0.874G 

0.3714 

31.238 

0.8401 

225 

0.8919 

0.3728 

30.449 

0.S33S 

230 

0.9099 

0.3744 

29.671 

0.8274 

240 

0.9484 

0.3775 

28.145 

0.8140 

250 

0.9907 

0.3807 

26.655 

0.7993 

260 

1.0385 

0.3845 

25.196 

0.7852 

270 

1.0905 

0.3879 

23.762 

0.7679 

275 

1.1198 

0.3900 

23.052 

0.7599 

28  Of 

1.1502 

0.3919 

22.347 

0.7506 

290 

1.2185 

0.3962 

20.947 

0.7321 

900 

1 . 2979 

0.4010 

19.557 

0.7126 

310 

1.3900 

0.4059 

18.170 

0.6905 

320 

1.5003 

0.4116 

16.782 

0.66G9 

325 

1 . 5625 

0.4144 

16.0S6 • 

0.6533 

330 

1.6315 

0.4173 

15.388 

0.6392 

340 

1.7933 

0.4237 

13.982 

0  60S6 

350 

1.9936 

0.4300 

12.558 

0.5721 

The  numbers  show 

n  above  r 

ipply  to  the  liquid  pi 

lase  along  the 

boiling  curve 

water. 


Values  of  o,  for  Some  Individual  Ions  in  Aqueous  Solutions 

o,  x  I0e 

Ion 

2.5 

3.0 

3.5 

4.0-4.S 

Rb\  Cs  +  ,  NHJ,  Tl  +  ,  Ag  + 

K  +  ,  CI-,  Br-,  1",  NOj 

OH",  F-,  HS",  BrO;,  IO;.  MnO; 

Na  +  ,  HCO;,  H2PO;,  HSO;,  Hg 2  + ,  SO;-,  SeO;",  CrO;- 
HPO;  ",  POJ- 

4.5 

5.0 

6 

'i® 

9 

II 

Pb  +  +  ,  CO;-,  SO;-,  MoO;- 

Sr++,  Ba++,  Ra++,  Cd++,  Hg++,  S”,WO;- 

Ll\  Ca++,  Cu++,  Zn++,  Sn+  +,  Mn++,  Fe++,  Ni+  +,  Co+  + 

Mg++,  Be+  + 

H  +  ,  Al3-1,  Cr3  +  ,  trivalent  rare  earths 

Th4\  Zr1  +  ,  Ce1  +  ,  Sn4  + 

SOURCE:  Adapted  from  I.  M.  Klotz,  Chemical  Thermodynamics.  Enjlewood  Cliffs,  N.J., 
Prtnticc-Hall.  1950.  p.  331. 
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DAT  A  X/3Q0.  , -1.1 .39,-8.06,-7.72,  64.2 6 L-7.C 
I .  0  3  ,  1 . 3  3  t  9 . 1  0  6  / 

R F 4r?(Slr':)°HMT>i,  PHM A  X  ,  F02 M  I N  ,  FO  2MC  X 


,-2.59,-< 


i\PH=  (  (  PH M  A  X-PHM  T  It ) /'■ 5  ) -t- 1 
MFC!  2=  (  (  f°?HA  x-F  1RH  i '  ■  )  / r' .  2 r>  )  + 1 
PH(  1  )  = PH MIN 
F02(  1  )  =  F G 2 Y  I N 


3  9 ,  i" 


DO  1  J=?,NPH 
3  PH  f  J  )  -r>H  (  J  -  1  )+■'•.  F 

DO  ?  J  =  2  »  N.FG2 
?  FOR  (  J  )  =F  02  f  J-l  )  +  .25 
NTf/iTA  L=NPH»MFO? 

Jt'  =  0 

DO  A  J=1,MPH 

-  np  ^  k,  -  1  ,  hr _ 

jo  jr  4-: 
x° ( jc ) -ph( j ) 

6 _ CPHTjni  ic  ___ _ 

?  TnfTf’iNUF 


N  C  ~  C 

no  7  j  - 1,  t  o  o'  i 

DO  o  K=1 »Mr02 
NC=NC+1 
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.  A  *  F  S  7  )  +  (  A 

L  u  ■ » 

i 

i  (  Y  (  ?  ) 

)  -  (  (  ?  .  AH  )  +X  (  9  )  4-  (  Q  .  5*f02  (  K  )  )  )  )  ) 

.  A !  1 

S  - 

i 

A 

!  i  4  A  SJ.-Xi.AJ 

_ _ _ _ _ _ - . . . . -  .  . . j 

AH2  S - ( AH+AHS ) -X { 3 ) 

ASU4=X  (  5  )  +  AS  +  (  2  -  F  (K  )  ) 

A  H  S  C  4  =  (  AH-*- A  SO  4)  -  X  (  6  ) 

AH?  s r 1 4  —  (  AH4-/\Hsn^  )  -v  (  1  s  ) _ 

QH2S04  =  1 0 ** ( A M2  S  04  - v { o)  ) 

X 1 0  (  K  C  )  =  A L  00 l r ( OH? 504 ) 

QS.fJ  (  A  S -  Y  (  7  )  J . . . 

X?  (  KG)  =  A L 0 G 10  C  OS) 

QHS=  1  O5"-*  (  AHS-  Y (  4  )  ) 

x? ( kg )  =  alt oio ( ohs  ) _ 

QH2S  =  1  ’**  (  AH?  S-y  (  o)  ) 

XI  (  KC)=Al.OGIO(OH2S) 

0504  =  10**  (  A  <:  '/, -v  (  /  )  )  . 

X 4 ( K C ) - A  I  HO 10 ( USO* ) 

QHS 0-4  =  1  *  *  (  AH  SO 4  -Y  (  5  1  ) 

X  ^  (  X  f  )  ~  A  1,-^Ql  0  (  OHS  ) _ 

Tl=  (  2  .  -  Y  (  1  )  -  (  OH  S+OHS04  +  4 . 4  (  QS  -*-05 04  )  )  ) 

T 2-  (  QOH  +  Q H5  +  0HS044-?  .  *(  9S+QS04  )  -Oh  ) 

QN  A=  (  P .  ?*  ( T  1  +T2  )  )  /  (  1  -*- P . ) _ 

X  1  5  {  K C ) ~ 00 A 
QK  =  R*CN!A 

y i 4 f or ) =0* _ 

AKS04=  (  (  \  LOG  1  0  (  OK, )  )  +v  (  ?  )  +  A  S04  }  -X  (  7  > 

A  NS  04  =  <  A  L  0  G 1  o  (  0 '  i  A  )  +Y  (  A  )  +  AS04  )  -  X  (  7  ) 

A  K  H  S  ( '  -4  =  (  \\  10  01  A  (  )  +  Y  (  ~  )  +  A  H  S04 )  -  y  (  F  )  _ _ _ 

A0H50 4=  ( Tal no  1 :  (  Q r ; A  )  +  Y  {  5_i  +  A H SO 4  j  -  X (  8  )  j 
0  X.  5  0  4  - 1  "  *  (  A  K  5  n  4  -  v  (  4  )  ) 

X  M  X  r  )  =  A L  00 1  i  (  04  5 n--f  ) _ 

QNSf  4=10*-' (  A  W4-Y{  0  )  } 

X7(  KC  )  =  -M_  AO 3  o  (  QMS 04  ) 

OKHSP4-3  C  A  XHS04  -  Y(  5  )  ) _ _ 

X  rt  (  K  C)  =  A  L  10 1  0  (  OK  MS  G-o  ) 

ONH  SP4  =1  G**  (  A  NH  S  n4 -Y  (Si) 

X  9  f  *  c  )  =  A  l  r;o  1  0  (  C  N't :  4  o 4  ) _ 

T  S 1  -  G  H  ?  S  +  0HS4-OS 

T  S  2=  O  sn  4  +  0  H  S  n  4+ A  K  S04+  QN  S  04  +  OK  HS  04  +  ON  H  S  Cl  4+QH2S  04 

TS=TS i+TS2 _ 

X  12 ( KC ) =TS 
XH2  S  =  OH?  S/T  S 

XHS  =  QHS/T  S _ 

XS=G'S/TS 

XSn4=0S04/TS 

X H  S  04-0  H S  Q4 / T  $ _ _ 

XKS04  =QKS04/TS 
XNS04=0NSn^/TS 

y  14 1 ;  4  r  4-0404^^/1  s  _ _ 

X 1 : tJ  S  04=  0  N  H  S  O  4  /  T  S 
XH2Sf)4-=QH2SU4/TS 

XRS-TS1/TS _ 

X  0  S  -  T  S  ?  /  T  5 
X  1 1  (  K  C  )  =  X  n  S 

r>  V4  S=-l  »  (  (  v  H  ?  S  X  M  1  )  )  -*-  (  XH4*X  (  1  2  )  )  +  (  xns*X  (  H)  )  ) 

X  i  8  (  K  C  )  = : '  5  ?  4  S 

OCl=  (  QHA  +0K.  )  -  (  OK  SO 4  +0NS04  +  T?  ) 

XI  6( KG) =OCL _ 


T  I  ?  = 0  .***<  Q M  A 4  Q K  +  0 C  L  +  DHS  +  OH 5 0 4  +  0 K S 0 4  +  0 N S  . : 4+4T  ( QS+QS04 )  ) 
X  1 7  f  K .  C  )  =  T  T  ? 

F  H  2 S  = ( AL 0010 ( 0H2S )  ) +  < ALOGin Z( 1 )  5  )-( ALPG10 (55.5)  ) 

rs= ( fhzs  +  (n .  i?  (k  )  n -  ( 7  m  +  Ai  non-  ( 7 ( 2  m*2 _ 

X  1 3  (  K  f  )  -  F  s 
KC  =  KG  +  1 
CONTI  KM  .IF 
3  CON? I NUF 


wp  nr(  , .  i  •-}  ) 


101  Ff  P.F  AT { / / 2  X  »  •  TEUPPPATUPF  KH20  KH2S  KHS  KSC4 

.  KHS 04  K  K S 04  »  { K  NAS 04 )  KHS 0 4  { NAHSD4 )  ' /?X , * 

K  J=r F  S2  A r IP  H 0 0  KH2SQ4  DEL  H2S-S 

• ; j  1 1  H S  —  S  1 1 i_L  SQ4-S  1  /  /  ) 

WP .  I  T  G  (  6 , 1  j2  )  {  X  (  !  )  ,  X  (  ?  )  ,  X  { 3  )  ,  X  (  4  )  ,  X  (  5  )  ,  X  {  A  )  ,  X  {  7  )  ,  X  {  8  )  ,  X  (  °  )  , 

.  x  ( i  ),v(n),y(i?),x(i?>) 


i  0  F('PM  AT  (  1  3 no  .  3  ) 

WP  I  T  E  (  f  ,  1  r  3  ) 

103  F  ( WM  ATC  //2X  «  1  ION STRENGTH  ACT  HOC  . ACTI  VI  TY COEF  . 

.OF  {  K  ,  C  L  )  (OH,  HS)  <  HS04,NA  ,  KS04  ,NA  SCJ4,  KHS04,  NAHS04  )  G/2^  1 
.504  S--  H  (H25, H2S04 ) ’// ) 

N  p I T  F  ( 6 ,  1 04  )  {  V (  1  )  , Y  m  1 Y (3)  , Y ( 4 )  ,  v ( n  , Y  f  6 }  , Y ( 7 )  , Y  M  )  ,  V  <  q )  ) _ 

104  F  0  k  4  a  t ( q  f 1 C . 4 ) 

VP.  I  T  E  (  6  ,  1  :  ?) 

1 0  5  FORM  A  TJ  /  /  2  X 1  HFNR  Y  CONSTANT  OF  H2  S . FIJGACI  TV . QF . WATER 

•HOLLAND  DATA  OF  H2S  *//) 

vn  to  f  n  n n  ( 7  ( i )  T  z  (  n  ,  z  ( 3  >  > 
m  FH;^Tnry,rio,oi  x,rOi,lV,Fir.2  ) 


K/NA  RATIO 


FUG AC  I T Y  OF  SULFUR  •//) 


WRITE (6, 107) 

137  FOR 7  AT ( / / 2 X ,  * 

_ V  K I  t  r  (  6  j  10 HJ  (  Ft  FS2  ) 

1 0 0  F 0 RMA T ( 1)  X , F 1 o . 3 ,  ] 0 X,F1C.4) 

WRITE (6,200  ) 

W”!  Tr  (4  t  ,?  n  )(>'?('),  VP  00  (  v  )  ,  X  1  (  ■  i  )  ,  X?  (  I- ■  )  ,  x  3(  N  )  t  X4(  N  )  ,  X5  (  N  )  ,  XM  r:  )  f 

.  X7c: )  , x*  (  • )  ,*9(0 ) ,  xu,  ( u ) ,  xi  1  ( n  ,  r:=i , NT  n  al  ) 

WRITE (6, 204) 

_ WR I TF ( 6 , ?02 ) 


WRITE  (  A  ,  2  03  )  (  XP(  M)  ,  X  F  02  (H)  ,  X  12  (  H)  ,X13(N)  ,X14(N)  , XI 5 ( N ) , X16 ( F  )  , 
. X 17 ( E  *  ,X I  3 ( N )  , N=1 , NTOTAL 5 


1 0  FORMAT  (  4F1D.2  ) 

3  0  p  F  0  P :  ‘  AT  (//.?*,  ’  PH'  ,  3X  ,  ' 17  02  »  ,  6X  ,  1  KHZ  S'  ,  7X  ,  *  Xh  S  ’  ,9X, • XS*  ,  BX,  1  *SU4»  ,  7  X 

•  1  XH$Q4>  t  6X,  1  XKS04 1  ,  5X  ♦  1  XNASO 4.1  ,  5 X,  1  XKHS04 1  ,  5X,  1  X N A H S 04  * . . 

4 X ,  1  XH2SG4 ' , 6X , 1 XO S  * //  ) 

201  F  0  Rjv  AT  ( r  5 . 1  ,  IX,  F  6.2,  IX,  F  10. 3,  1*  ,  F  1 0  .  3,  1  X, FI  0.3, IX, FI  0.3,  IX, FI  0 . 3  , 

,  i  v  n  i  n  3 ,  i  x .  r  i  c  n ,  i  * ,  n  ~  ♦  3 ,  i  t  ,  f  i  o .  m  x  ,  n  o .  n  i  k  ,  -  i  o .  n 


fi  o 


FORMAT  (  //2X  ,  *  P H  '  ,  4X  ,  ' F ( : 2  *  ,  7 X ,  'IS*  ,  1 1  X  ,  '  F  S  '  ,  9 X  ,  1  MK  ’  3 X  ,  1  .MN  A  *  ,  9 X  ,  '  M C L 


.  •  , 14 y  ,  » ] 2*  ,ox,  ' 9 S3 4  S * // ) 

203 . FCI P.H AT  (  F 5 . 1  , 1 X ,  F 6 .  2 , 2 X  «  E  1 0 . 3 .2*  ,  FI C  .  3  X  ,.F  19 . 7  ,  °X  ,  "1  f  F  ,.2i.t 

. F 10.3,2X,Fln.3,2X,F10*3) 

2  n  F '  1 P  I A  T  (  1  1  ’  ) 


